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SUMMARY 


The  support  provided  by  the  grant  AFOSR-90-0210  has  resulted  in  the  publication 
of  four  papers  and  three  conference  proceedings.  Our  research  covered  three 
topics . 

First  we  studied  shock  waves  generated  by  dielectric  breakdown  of  nano-  and 
picosecond  pulses  at  1.06  fm  wavelength  in  water  [P-1].  We  described  a  non- 
invasive  optical  technique  to  measure  the  propagation  velocity  of  shock  waves  in 
water  and  showed  that  during  spherical  expansion  in  water  the  pressure  of  a 
shock  wave  scales  as  the  inverse  of  the  square  of  the  propagation  distance  and 
drops  to  acoustic  values  within  a  few  hundred  microns. 

Second  we  investigated  the  pressure  transients  generated  during  the  ablation  of 
sheet  polyimide  by  nanosecond  excimer-laser  pulses  [P-2,  C-1,  C-2].  We  observed 
the  formation  of  multi-kbar  shock  waves.  The  pressures  generated  during 
ablation  were  shown  to  result  from  the  violent  expansion  of  the  rapidly  heated 
gaseous  ablation  products.  Initially,  this  work  used  an  optical  system  to 
measure  shock  wave  velocities.  The  use  of  this  system  was,  however,  technically 
quite  demanding.  We  subsequently  investigated  the  use  of  x-cut  quartz  gauges  to 
measure  shock  pressures.  The  initial  measurements  using  these  gauges  showed 
that  their  sensitivity  was  insufficient  to  measure  acoustic  pressures  at  or 
below  ablation  threshold. 


We  then  developed  a  system  based  on  PVDF  (polyvinylidene  fluoride)  thin-film 
transducers  which  allowed  measurement  of  the  pressure  transients  generated 
during  excimer  laser  ablation  of  polyimide  for  fluences  between  30  J/m^  and 
1  MJ/m^.  Measurements  were  made  for  four  wavelengths:  193,  248,  308,  and 
351  nm.  Three  distinct  fluence  ranges,  within  which  separate  physical 
mechanisms  govern  pressure  generation,  were  found.  In  the  low  fluence  regime, 
the  pressure  is  believed  to  be  due  to  subsurface  thermal  decomposition  for  the 
two  longer  wavelengths  and  to  ablative  photo-decompositions  at  the  two  shorter 
wavelengths.  At  higher  fluences  the  pressures  result  from  thermal  expansion  of 
ablation  products  and,  at  sufficiently  high  fluences,  from  plasma  formation. 

The  measurement  technique,  demonstrated  for  polyimide,  also  has  potential 
application  to  biological  materials  such  as  bone,  skin,  and  cornea. 


Third,  in  addition,  work  initiated  prior  to  the  AFOSR  grant  by  ADZ  at  the 
University  of  Berne,  Switzerland,  was  completed  [P-3,  P-4,  C-3,  C-4,  C-5].  That 
research  was  concerned  with  the  ablation  of  tissue  and  tissue-like  targets  by 
pulsed  laser  radiation  in  the  infrared.  A  model  that  explains  the  ablation 
process  as  a  combination  of  evaporation  and  ejection  of  liquefied  material  was 
developed . 


Publications  resulting  from  this  grant  are  attached  and  describe  the  research  in 
detail . 
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Abstract  Shock  waves  generated  by  a  laser-induced  plasma  were  investigated  using  a  pump- 
and-probe  technique.  Both  7-ns  and  40-ps  laser  pulses  at  1.06  pm  were  employed  to  initiate 
breakdown  in  water.  Two  He-Ne  laser  beams  were  used  as  a  velocity  probe,  allowing  the 
accurate  measurement  of  the  shock  velocity  around  the  plasma.  The  maximum  shock  pressure 
was  determined  from  the  measured  shock  velocities,  the  jump  condition  and  the  equation  of 
state  for  water.  The  conservation  of  the  total  momentum  of  the  shock  front  was  used  to  derive 
expressions  for  the  shock  velocity,  particle  velocity  and  shock  pressure  vs.  the  distance  (r) 
from  the  center  of  the  plasma.  For  a  shock  wave  of  spherical  symmetry,  the  shock  pressure  is 
proportional  to  1/r^.  Our  work  shows  that  the  expanding  plasma  initially  induces  a  shock  wave; 
the  shock  wave  dissipates  rapidly  becoming  an  acoustic  wave  within  300-500  pm. 

PACS:  42.80.-f,  87,  43.35.Cg 


The  interaction  of  laser  radiation  with  matter  can  lead 
to  the  generation  of  pressure  waves.  Depending  upon 
the  type  of  interaction,  the  generated  waves  can  be  ci¬ 
ther  acoustic  waves,  low-pressure  waves  propagating  with 
the  speed  of  sound  or  shock  waves,  high-pressure  waves 
propagating  at  supiersonic  speed  [1].  The  latter  appear 
when  the  absorption  of  laser  radiation  is  followed  by  a 
rapid  phase  change  of  the  medium  such  as  evaporation 
or  formation  of  plasma  [2,3]. 

Our  interest  in  the  study  of  shock  waves  stems  from 
an  increasing  awareness  that  the  pressure  waves  gener¬ 
ated  during  the  medical  use  of  lasers  may  have  impor¬ 
tant  effects  on  tissue.  Laser-induced  optical  breakdown 
is  currently  used  in  two  medical  applications;  kidney 
stone  fragmentation  in  urology  [4]  and  photodisruption 
in  ophthalmology  [5].  In  the  former  application  the  de¬ 
sired  fragmentation  of  the  stone  is  most  likely  caused 
by  shock  waves.  In  photodisruption  the  desired  effect  is 
probably  caused  by  the  plasma  itself;  the  ancillary  shock 
waves,  however,  have  the  potential  of  causing  damage  to 
cells  and  subcellular  structures  [6,7]. 

The  high  peak  intensity  of  short  duration  laser  pulses 
inside  a  transparent  medium  leads  to  optical  breakdown. 
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In  the  picosecond  to  nanosecond  time  regime,  the  break¬ 
down  is  initiated  via  a  multiphoton  avalanche  mechanism 
[8].  The  process  begins  at  the  focus  of  a  laser  beam,  where 
the  laser  irradianoe  is  the  highest  In  water,  optical  break¬ 
down  is  followed  by  the  formation  of  plasma.  The  high 
absorbance  of  the  plasma  facilitates  the  deposition  of 
laser  energy  in  the  focal  zone  even  in  transparent  media. 
Once  generated,  the  hot  plasma  expands  adiabatically 
with  supersonic  velocity  [9,  10],  creating  a  high  pressure 
front.  As  the  gases  in  the  cavity  cool,  the  expansion  veloc¬ 
ity  of  the  plasma  rapidly  decreases  and  the  high  pressure 
front  (shock  wave)  separates  from  the  cavity  or  bubble. 
The  expansion  of  the  bubble  continues  until  all  the  initial 
kinetic  energy  is  transformed  into  the  potential  energy  of 
the  cavitation.  Subsequently,  the  bubble  implodes  adia¬ 
batically  under  the  external  pressure  of  the  liquid.  The 
temperature  and  pressure  inside  the  bubble  rise  again 
leading  to  a  second  explosion  and  the  formation  of  a 
second  shock  wave  [1 1]. 

A  primary  characteristic  of  a  shock  wave  is  the  fast 
rise  time  of  the  pressure  wave,  which  for  all  practical 
purposes  amounts  to  a  discontinuity  in  pressure,  density, 
particle  velocity  and  internal  energy  [12].  In  water,  the 
rise  time  of  a  shock  wave  is  of  the  order  of  a  picosecond 
1 1 3],  corresponding  to  a  shock  front  thickness  of  a  few 
nanometers.  The  large  pressure  gradient  and  the  high 
translational  velocity  of  the  molecules  within  the  shock 
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front  account  for  the  unique  interactions  of  the  shock 
waves  with  matter  [14,  15]. 

Since  the  first  measurements  of  Bell  and  Landt  [3],  a 
number  of  different  techniques  have  been  developed  for 
the  study  of  shock  waves.  Pressure  measurements  using 
transducers  [16]  have  been  augmented  by  optical  meth¬ 
ods  such  as  Schlieren  photography  [17],  fast  photography 
[18]  and  optical  pump-and-probe  techniques  [19-21].  The 
advantage  of  the  optical  techniques  is  that  they  are  non- 
invasive.  In  addition,  the  spatial  resolution  is  potentially 
high,  of  the  order  of  micrometers,  and  the  time  resolution 
is  only  Umited  by  the  bandwidth  of  the  detector,  currently 
available  in  the  hundreds  of  GHz  range,  corresponding 
to  a  time  resolution  of  a  few  picoseconds. 

In  previous  pump-and-probe  experiments  [10,  20]  the 
initial  shock  velocity  was  calculated  from  the  simulta¬ 
neous  measurement  of  the  radial  distance  traveled  by 
the  shock  wave  and  the  time  of  arrival.  We  have  em¬ 
ployed  two  probe  beams,  separated  by  a  few  tens  of 
micrometers,  to  measure  the  shock  velocity  [22].  This 
experimental  arrangement  is  insensitive  to  the  temporal 
and  spatial  fluctuations  in  the  launching  of  the  shock 
wave.  From  the  traversal  time  of  the  shock  front  across 
the  two  probe  beams  and  the  distance  between  them 
we  derived  the  shock  velocity.  Using  this  high  resolution 
apparatus  we  have  mapped  the  shock  velocity  near  the 
plasma,  and  used  the  measured  shock  velocities  together 
with  the  jump  conditions  and  the  equation  of  state  for 
water  to  determine  the  maximum  shock  pressures  from 
first  principles.  We  have  applied  the  conservation  of  total 
momentum  of  the  shock  front  to  derive  an  expression  for 
the  maximum  shock  pressure  as  a  function  of  the  dis¬ 
tance  (r)  from  the  center  of  the  plasma.  We  have  found 
that  the  maximum  pressure  is  proportional  to  1/r^. 

An  important  aspect  of  shock  waves  as  far  as  medical 
applications  are  concerned,  is  the  range  over  which  they 
can  propagate  before  they  dissipate  enough  energy  to 
become  acoustic  waves.  In  this  work,  we  have  determined 
that  the  range  of  plasma-induced  shock  waves  is  limited 
to  a  few  hundred  micrometers.  In  comparison,  acoustic 
waves  have  a  range  of  several  tens  of  millimeters. 


1.  Experimental 

The  experimental  apparatus  is  shown  in  Fig.  1.  A 
Q-switched  (Schwartz  Electro-Optics)  and  a  mode-locked 
(Quantel  model  501)  Nd  :  YAG  laser  were  used  to  produce 
plasma  in  the  sample.  Both  lasers  operated  at  1.06  pm  anu 
produced  pulses  of  approximately  7  ns  and  40  ps,  respec¬ 
tively.  The  energy  of  the  lasers  was  monitored  with  a  Sci- 
entech  (model  365)  or  a  Molectron  (model  J3-09)  energy 
meter.  The  pump  pulse  was  focused  by  a  12-mm  focal 
length  aspheric  lens.  The  beam  waist  in  air  was  estimated 
using  burn  paper.  It  was  of  the  order  of  50  pm,  yielding 
an  irradiancc  in  the  range  of  10'®  to  10"  W/cm^  for  the 
nanosecond  and  10"  W/cm"  for  the  picosecond  pulses, 
respectively.  The  beam  waist  in  water,  however,  would 
be  bigger.  The  sample,  distilled  and  deionized  water,  was 
cjrculatecl  and  liltered  in  order  to  remove  bubbles  ami 
any  other  impurities  generated  during  the  breakdown 


Fig.  i.  The  expenmenlal  apparatus.  Two  He-Ne  laser  beams  are 
used  to  determine  the  velocity  of  the  shock  wave 


Two  He-Ne  laser  beams  separated  by  30-40  pm  were 
used  to  measure  the  shock  velocity.  The  system  of  the 
directing  mirrors,  focusing  and  collecting  optics  for  the 
probe  beams  was  placed  on  an  —  Z  translation  stage, 
X  being  the  axis  of  the  pump  beam  (axial)  and  2,  the 
axis  perpendicular  to  it  (radial).  A  dove  prism  in  the  path 
of  the  two  He-Ne  beams  was  used  to  make  the  plane 
of  the  two  beams  either  horizontal  or  vertical,  depending 
on  whether  the  measurements  were  carried  out  along  the 
axial  or  radial  direction.  The  two  He-Ne  beams  exiting 
from  the  cell  were  recombined  and  directed  onto  a  single 
Si  photodiode  (EG  &  G  model  FND  100).  The  active  area 
of  the  detector  was  5  mm^  and  the  rise  time  <  1  ns.  The 
signal  was  recorded  by  a  transient  digitizer  (Tektronix 
7912  AD,  750  MHz  bandwidth)  and  transferred  to  an 
IBM-AT  PC  computer  for  processing. 

The  ‘•hock  wave  induces  a  discontinuity  in  the  index  of 
refraction  of  the  medium  due  to  the  density  increase  in  the 
shock  front.  This  results  in  the  deflection  of  the  He-Ne 
laser  beam,  and  the  intensity  of  the  light  transmitted  to 
the  detector  is  decreased.  As  the  shock  wave  propagates 
through  the  two  probe  beams,  it  generates  two  separate 
pulses.  On  the  other  hand,  the  passage  of  the  cavitation 
bubble  through  the  probe  beams  produces  a  stepped 
signal,  resulting  from  the  interruption  of  one  and  then 
both  beams.  The  shock  velocity  was  determined  from  the 
time  interval  between  the  two  pulses  and  the  separation 
distance  of  the  two  probe  beams.  In  addition,  the  time 
interval  between  the  two  steps  was  used  to  determine 
the  expansion  velocity  of  the  plasma  or  the  cavitation 
bubble. 

Figure  2  shows  a  series  of  photodiode  signals  taken  at 
different  distances  from  the  center  of  the  plasma;  (a)  at 
the  center  of  the  plasma  where  both  beams  were  blocked, 
(b)  and  (c)  where  the  shock  wave  is  initially  launched  by 
the  expanding  plasma,  (d)  in  the  proximity  of  the  break¬ 
down  region  after  the  shock  wave  had  separated  from 
the  bubble,  while  (e)  shows  the  photodiode  signal  in  the 
l;ir  tielii.  The  wavelbmis  presented  in  big.  2h  and  d  cor¬ 
respond  l(»  the  stages  ol  the  pla.sina  e\'ohiti<’n  depicted 
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Fig.l  The  photodiode  signal  taken  at  different  positions:  (a)  the 
center  of  the  plasma,  (b)  the  edge  of  t’  e  plasma,  (c)  the  initial 
separation  of  the  shock  wave  from  the  plasma,  (d)  in  the  near 
field  and  (e)  in  the  far  field.  Rapid  events,  like  the  deflection  of 
the  beams  by  the  spatially-thin  shock  wave  front,  result  in  two 
separate  pulses.  Long  term  events,  like  the  passage  of  a  cavitation 
bubble,  produce  a  stepped  signal,  resulting  from  the  interruption 
of  one  and  then  both,  beams.  The  time  between  the  two  pulses  is 
related  to  the  shock  velocity,  while  the  time  interval  between  the 
two  steps  that  follows  is  related  to  the  velocity  of  the  expansion  of 
the  cavitation  or  bubble.  The  increased  signal  at  the  be^nning  of 
the  trace  is  caused  by  the  emitted  light  from  the  plasma 


in  Fig.  3b  and  c.  The  cavitation  bubble  can  undergo  sev¬ 
eral  cycles  of  nearly  adiabatic  expansion  and  contraction. 
Figured  illustrates  this  for  the  case  when  three  cavitation 
bubbles  are  formed.  Note  that  the  second  and  third  bub¬ 
bles  produce  half  as  large  a  signal  as  the  first,  indicating 
that  they  have  not  propagated  to  the  second  probe  beam. 

The  shape  of  the  photodiode  signal  was  an  almost 
symmetric  pulse  of  10  ns  (FWHM).  The  rise  time  cor¬ 
related  to  the  traversal  time  of  the  shock  front  through 
the  waist  of  the  probe  beam.  The  decay  of  the  signal, 
however,  could  be  altered  by  changing  the  diameter  of 
the  aperture  in  front  of  the  photodiode  We  found  that 
the  best  time  resolution  could  be  obtained  by  keeping  the 
aperture  open.  In  this  case,  the  diameter  of  the  collecting 
lens  is  the  effective  aperture.  Under  these  conditions,  the 
photodiode  signal  does  not  represent  the  course  of  the 
pressure  wave  in  our  experimental  arrangement.  Decay 
times  of  shock  waves,  reported  in  the  literature,  vary 
from  10  to  40ns  [21,  23].  A  possible  explanation  for  the 
short  duration  of  the  signal  in  our  experiments  is  that  the 
deflection  of  the  beam  is  highest  when  the  incident  beam 
is  at  grazing  angle  to  the  spherical  shock  front.  Therefore, 
only  the  leading  edge  of  the  shock  wave  causes  a  large 
deflection  that  can  be  detected. 

The  photodiode  signal  decreases  by  a  factor  of  three 
between  the  near  and  far  fields.  The  signal  intensity  de¬ 
pends  on  the  pressure  of  the  shock  wave.  In  this  work, 
wc  did  not  try  to  calculate  the  shock  pressure  from  the 
intensity  of  the  sijinal.  The  change  in  the  index  of  re- 


Fig.  X  A  schematic  representation  of  the  shock  wave  generation 
during  optical  breakdown:  (a)  plasma  initiation,  (b)  plasma  expan¬ 
sion,  (c)  separation  of  the  shock  wave  from  the  plasma,  (d)  collapse 
of  the  bubble 


Fig.  4.  Oscilloscope  trace  of  the  probe  signal  caused  by  a  series  of 
three  expansions  and  contractions  of  a  cavity.  Each  collapse  of  the 
bubble  is  followed  by  the  generation  of  a  shock  wave  as  shown  by 
the  arrows.  The  second  and  third  bubbles  produce  half  as  'srge  a 
signal  as  the  first  because  they  have  not  propagated  to  the  second 
probe  beam 

fraction,  caused  by  the  pressure,  is  proportional  to  the 
pressure  rise  only  up  to  2kbar,  for  higher  pressure  val¬ 
ues  it  increases  less  than  linearly  [19,  24].  These  types  of 
measurements,  however,  require  an  accurate  calibration 
of  the  detector,  which  is  not  necessary  when  the  measured 
quantity  is  the  time  between  the  two  pulses. 

The  boundaries  of  the  plasma  were  determined  us¬ 
ing  the  probe  beams.  The  probe  beams  were  translated 
through  the  plasma  to  the  position  where  both  beams 
were  blocked.  The  positions  for  total  blockage  of  the 
probe  beams  in  cither  direction  gave  the  plasma  ex¬ 
tremes.  The  mean  values  of  the  two  positions  in  the  axial 
and  radial  directions  were  taken  as  the  coordinates  of 
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the  geometrical  center  of  the  plasma.  It  should  be  noted 
that  the  center  of  the  plasma,  as  defined,  may  not  coin¬ 
cide  with  the  focal  point  of  the  pump  beam  [25]  or  the 
center  of  curvature  of  the  expanding  shock  wave.  This, 
however,  is  more  likely  to  be  significant  in  the  axial  than 
in  the  radial  direction.  We  will  return  to  this  point  in  the 
discussion. 

The  pressure- wave  velocity  was  measured  radially  and 
axially  as  far  as  1500  pm  away  from  the  center  of  the 
plasma,  on  both  sides  of  it.  The  range  was  chosen  so  that 
there  was  a  sufficient  number  of  measurements  in  the 
far  field  where  the  pressure  wave  propagates  with  sound 
velocity.  In  fact,  the  constancy  of  the  measured  velocity 
in  the  far  field  was  taken  to  imply  that  the  sound  velocity 
had  been  reached.  The  distance  between  the  two  He-Ne 
probe  beams  was  determined  from  the  traversal  time  in 
the  far  field  and  the  sound  velocity,  1.483  km/s  in  distilled 
water  at  20°  C.  The  distance  between  the  two  beams 
was  typically  30-40  pm.  Close  to  the  plasma,  the  shock 
velocity  was  measured  in  steps  of  20  pm.  In  the  far  field, 
where  changes  in  the  velocity  were  small,  the  velocity 
of  the  pressure  wave  was  measured  at  intervals  of  50  or 
100  pm.  There  were  considerable  fluctuations  of  the  shock 
velocity  in  the  proximity  of  the  plasma,  the  standard 
deviation  of  the  measured  velocities  ranging  from  60  m/s 
to  140  m/s.  The  velocity  fluctuations  were  attributed  to 
the  statistical  nature  of  the  shock  wave  generation  and  the 
variation  of  the  pulse  energy.  The  velocity  fluctuations 
decreased  gradually  with  increasing  distance  from  the 
plasma,  correlating  with  the  dissipation  of  the  shock  wave 
and  the  appearance  of  the  acoustic  wave.  The  standard 
deviation  of  the  measured  velocity  in  the  far  field  provides 
a  measure  of  the  velocity  resolution  of  the  apparatus. 
The  standard  deviation  of  the  measured  sound  velocity 
ranged  between  15  m/s  and  20m/s,  corresponding  to  a 
resolution  of  the  experimental  apparatus  of  ~  30-^ m/s 
or  2-3%. 

The  plasma  duration  produced  by  the  nanosecond 
pulses  was  of  the  order  of  1 5  ns.  The  emission  spectrum 
of  the  plasma  was  collected  with  an  optical  fiber  coupled 
to  1 /4-meter  flat-field  spectrograph  and  measured  using 
a  gated  (30  ns)  optical  multichannel  analyzer  (EG&G 
OMA  III).  The  emission  spectrum  showed  a  broad  struc¬ 
tureless  band  with  the  maximum  intensity  at  ~  370  nm. 
Using  Wien's  law  (An,T  =  0.29  K  ■  cm),  we  have  calcu¬ 
lated  a  plasma  temperature  of  ~  8000  K  which  is  of  the 
same  order  of  magnitude  as  the  plasma  temperatures 
reported  in  the  literature  [9],  The  emission  spectrum, 
however,  deviated  from  that  of  black-body  radiation  in 
the  long  wavelength  region,  possibly  because  the  gate 
was  too  long  to  allow  temporal  resolution  of  the  chang¬ 
ing  spectral  emission  and  the  observed  spectrum  was  the 
superposition  of  the  spectra  emitted  during  the  evolution 
of  the  plasma. 

The  plasma  size  depended  on  the  energy  and  duration 
of  the  laser  pulse.  Nanosecond  pulses  with  energy  of  2 
3  mJ  generated  plasmas  that  were  spherical  or  ellipsoidal 
with  a  ratio  of  major  (axial)  to  minor  (radial)  axes  varying 
from  1  :1  to  1.5:1,  the  minor  axis  being  approximately 
80  lOOpm  When  the  cnergv  increased  to  11  15m.I.  the 
r.itio  of  the  axes  ranged  from  2:1  to  2  I.  .'in<l  the  minor 


axis  increased  to  120  150  pm.  The  breakdown  threshold 
was  1  -2mJ.  For  picosecond  pulses,  3  mi,  the  ratio  of  the 
two  axes  was  3:1,  the  minor  axis  being  approximately 
50  pm.  In  both  cases  the  plasma  dimensions  varied  from 
pulse-to-pulse  with  the  axial  dimension  varying  more 
than  the  radial. 


2.  Results 

Figure  5  shows  the  results  of  a  typical  measurement  of  the 
shock  velocity  produced  by  nanosecond  laser  pulses  and 
measured  axially  from  the  center  of  the  plasma.  Simi¬ 
lar  results  were  obtained  in  the  radial  direction.  Each 
experimental  point  is  the  average  of  ten  measurements 
and  the  error  bars  are  the  standard  deviations.  The  mea¬ 
sured  shock  velocity  varies  from  ~  2.2  km/s  at  a  distance 
of  ~  120  pm  from  the  center  of  the  plasma  to  that  of 
the  sound  velocity  at  500  pm  away  from  the  center  of 
the  plasma.  Our  measurements  in  the  proximity  of  the 
plasma  agree  with  the  velocity  measured  by  Fujimoto 
etal.  [20]. 

The  shock  wave  pressure  is  related  to  the  shock  veloc¬ 
ity  through  Newton’s  second  law  across  the  shock  wave 
discontinuity. 

P-Po=  U,Up(?o,  _  (1) 

where  P  and  Pq  are  the  shock  and  the  hydrostatic  pres¬ 
sure,  respectively,  Ug  the  shock  velocity,  Up  the  particle 
velocity  and  go  the  density  of  the  water  before  compres¬ 
sion.  Po  is  much  smaller  than  the  shock  pressure  and  can 
be  neglected.  Equation  1  and  the  equations  that  describe 
the  conservation  of  mass  and  energy  across  the  shock 
wave  discontinuity  are  known  as  the  jump  conditions. 

The  shock  and  particle  velocities  are  related  through 
the  equation  of  state.  To  a  first  approximation,  it  can  be 
written  as 

Us  =  A  +  Bup.  (2) 


AXIAL  DISTANCE  FROM  THE  CENTER  OF  THE  PLASMA  (mm) 


Fig.  5.  Shock  velocity  measured  axially  vs.  distance  from  the  center 
of  the  plasma.  The  shock  wave  was  generated  by  a  2.6.mJ  pulse  from 
a  (t-switched  Nd  :  YACj  laser.  The  arrow  indicates  the  direction  of 
the  exciting  la.scr  pulse.  Irach  experimental  point  is  the  average  of 
ten  measurements.  The  large  diflercnce  in  the  values  of  the  shock 
velocity  al  the  proximity  of  the  plasma  is  to  some  extent  accounted 
for  by  the  dillerent  nninnnnn  distances  lioin  the  center  plasma;  fi'i 
ex.imple  l2T|mi  al  the  proximal  side.  l.‘'0(im  al  the  di'-lal  side  1  he 
emu  h.irs  lepieM'ii!  llie  si.uid.iui  dexi.ition 
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AXIAL  DISTANCE  rSOII  THE  CENTER  OF  THE  PLASMA  (mm) 

Fig.  6.  Shock  peak  presssure  vs.  axial  distance  from  the  center  of 
the  plasma.  Each  point  is  the  average  of  ten  values  of  pressure 
determined  separately  from  the  shock  velocities  shown  in  Fig.S. 
The  error  bars  represent  the  pressure  standard  deviatioa  The  arrow 
indicates  the  direction  of  the  exciting  laser  pulse 


For  water,  the  values  of  A  and  B  are  1.48  km/s  and 
2.07,  respectively  [26].  t/,  and  Up  are  given  in  km/s.  Note 
that  the  value  of  A  is  the  sound  velocity.  These  values 
of  A  and  B  are  valid  for  up  to  20kbar.  A  higher  order 
approximation  is  given  in  [27].  From  Eqs.  (1)  and  (2) 
the  shock  pressure  is  calculate^  as  shown  in  Fig.  6.  The 
calculated  pressures  range  from  8kbar  (8  x  10*  Pa)  at  a 
distance  of  ~  100  pm  from  the  center  of  the  plasma  to  a 
few  hundred  bar  at  a  distance  of  500  pm. 

The  propagation  of  the  shock  wave  involves  the  in¬ 
elastic  collisions  of  the  water  molecules  within  the  shock 
front.  For  a  shock  wave  that  has  separated  from  the  ex¬ 
panding  bubble,  the  total  momentum  of  the  shock  front 
is  conserved.  For  a  spherically  symmetric  shock  wave,  the 
conversation  of  the  total  momentum  at  a  distance  r  from 
the  center  of  the  plasma  can  be  written  as; 

47rr^(ysdtpwp  =  const,  (3) 

where  At  is  the  rise  time  and  q  the  density  of  the  medium 
in  the  shock  front  The  density  within  the  shock  front 
changes  from  qq  to  that  of  the  compressed  state.  The 
density  of  water,  even  at  lOkbar,  increases  by  less  than 
15%  [27]  and  can  be  considered  constant  At  is  the  rise 
time  of  the  shock  front  and  is  constant  in  the  range  of  0- 
lOOkbar  [13].  Substituting  Up  from  the  equation  of  state 
of  water  (2)  into  (3),  we  derive 

U,{U,  -A)=~,  (4) 

where  C  includes  the  constants  B,  g  and  At.  Equation 
(4)  can  be  solved  for  Us  to  derive  the  expression  of  the 
shock  velocity  vs.  distance  from  the  center  of  the  plasma 

Equation  (5)  is  the  mathematical  expression  used  to  fit 
the  measured  velocities  of  the  shock  waves  generated  by 
nanosecond  and  picosecond  pulses  in  both  radial  and 
axial  directions.  The  constant  C  was  calculated  by  a 
weighted  linear  regression  using  the  measured  values  of 
the  velocity  vs.  distance.  For  the  purpo.sc  of  tilting.  (4) 
can  be  written  as  Y  =  CX  where  >'  -  liXl\  -1| 


and  X  =  1/r^.  The  value  of  C  is  not  sensitive  to  the 
measured  velocities  in  the  far  field.  In  fact,  the  difference 
in  the  value  of  the  constant  obtained  by  using  the  full 
range  of  measurements  (±1500  pm)  instead  of  a  reduced 
range  (±450  pm)  is  less  than  10%.  The  justification  for 
using  a  truncated  range  is  that  the  shock  wave  range  is 
of  the  order  of  500  pm. 

There  are  two  sources  of  error  in  the  measurement 
of  the  shock  velocity.  First,  the  calculated  shock  velocity 
is  a  statistical  average  of  many  separate  measurements, 
rather  than  the  velocity  of  a  single  generated  shock  wave. 
Second,  the  velocity  measured  is  the  average  velocity 
between  the  two  probe  beams  which  can  be  written  as: 

'HL/2) 

mr))  =  (l/L)  J  Us{r')dr',  (6) 

r-(L/2) 

where  L  is  the  distance  between  the  two  probes  and 
Us(r')  is  the  shock  velocity,  given  by  (5),  The  integral  can 
be  evaluated  but  the  derived  expression  is  complicated 
and  not  easily  amenable  to  a  standard  fitting  procedure. 
In  view  of  the  large  fluctuations  of  the  measured  shock 
velocities,  the  use  of  a  deconvolution  procedure  would 
not  add  to  the  accuracy  of  the  present  analysis.  In  fact, 
the  numerical  evaluation  of  integral  6  over  a  distance  of 
40  pm  (the  distance  between  the  two  probes)  gave  a  value 
that  was  within  the  experimental  standard  deviation  of 
the  measured  shock  velocity. 

Substituting  (5)  into  (2)  (equation  of  state),  we  derive 
an  expression  for  the  particle  velocity : 


A^  C 

“T  +  T 


(7) 


We  can  use  (7)  to  extrapolate  the  particle  velocity  at  the 
onset  of  the  shock  wave.  The  shock  wave  separates  from 
the  cavitation  when  the  particle  velocity,  that  is  the  ve¬ 
locity  of  the  displacement  of  the  matter  behind  the  shock 
front,  exceeds  the  velocity  of  the  bubble  expansion.  The 


Table  1.  Comparison  of  the  bubble  expansion  velocity  and  particle 
velocity  at  the  onset  of  the  shock  wave 


Energy 

[mJ] 

Distance 

lum] 
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velocity 

[km/s] 

Particle 

velocity 

[km/s] 

Pressure 

[kbar] 

Nanosecond-generated  shock  wave 
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-  60 

0.26 

0.52 

13.1 
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0.52 
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+  110 

0.49 

0.37 

8.4 

14  7  (R) 

-120 

0.50 

0.57 

15.3 

+  120 

0.49 

0.57 

2.6  (A) 

-no 

0.33 

0.32 

6.7 

+  102 

0.38 

0.42 

9.9 

11.3  (A) 

-203 

0.39 

0.23 

4.6 

+203 

0.38 

0.24 

6.1 

Picosecond-generated  shock  wave 

3.3  (R) 

-  50 

0.24 

0.48 

11.9 

1  50 

o.;9 

0.48 
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expansion  of  the  bubble  can  be  measured.  The  lime  in¬ 
terval  between  the  two  steps,  in  Fig.  2b,  corresponds  to 
the  traversal  time  of  the  expanding  plasma  between  the 
two  probe  beams.  From  the  traversal  time  and  the  dis¬ 
tance  between  the  two  beams,  we  calculated  the  velocity 
of  expansion.  Table  1  shows  the  measured  values  of  an 
average  velocity  of  the  plasma  expansion  and  the  cal¬ 
culated  particle  velocities.  In  most  cases  the  agreement 
is  very  good.  The  process  of  separation  of  the  shock 
wave  from  the  cavitation  described  here  is  limited  by  the 
temporal  resolution  of  our  experimental  apparatus.  The 
actual  process  may  be  much  more  complex. 

2.1  Radial  Shock  Velocity 

Figures  7  and  8  show  the  measured  radial  velocities  of 
the  shock  waves  generated  by  nanosecond  and  picosec¬ 
ond  pulses,  respectively.  The  lines  represent  the  theoret¬ 
ical  fit  of  (5)  to  the  data  of  both  branches.  Equation  (5) 
can  also  be  fitted  separately  to  each  of  the  branches.  The 
difference,  however,  is  within  the  experimental  error.  The 
theoretical  curves  agree  reasonably  well  with  the  mea¬ 
sured  shock  velocities,  considering  the  simplicity  of  the 
model  used.  The  same  expression  for  the  shock  veloc- 


AXIAL  DISTANCE  FROM  THE  CENTER  OP  THE  PLASMA  (mm) 

Fig.  9.  Shock  velocities  vs.  axial  distance.  The  shock  waves  were 
generated  by  2.6  or  11.3mJ  pulses  from  a  g-switched  laser.  The 
lines  are  the  theoretical  fit  as  given  by  (5).  Both  branches  of  the 
theoretical  curve  begin  at  the  same  distance  from  the  center  of  the 
plasma.  Each  branch  was  fitted  separately.  The  arrow  indicates  the 
direction  of  the  exciting  laser  pulse 


ity  vs.  distance  applies  to  the  shock  waves  generated  by 
picosecond  and  nanosecond  pulses.  In  this  respect  the 
pulse  duration,  at  least  in  the  range  10~“-I0“*s,  does 
not  affect  the  scaling  of  the  shock  velocity.  The  constant 
in  the  expression,  however,  does  depend  on  the  laser 
parameters. 


RADtAL  DtSTANCE  FROM  THE  CENTER  OF  THE  PLASMA  (mm) 

Fig.  7.  Shock  velocities  vs.  radial  distance.  The  shock  waves  were 
generated  by  3.6  or  14.7  mJ  pulses  from  a  ^-switched  laser.  The 
lines  are  the  theoretical  fit  as  given  by  (5).  Both  branches  were 
fitted  simultaneously 


radial  mSTANTE  FROM  THE  CENTER  OF  THE  PLASMA  (mm) 

l-ig.  8.  .Shock  velocities  v.s.  radial  distance,  rhe  shock  veaces  were 
cenerak'd  h)  a  ni.I  pulse  from  a  mode-locked  l.isei  [he  line 
IS  the  IhcoretK'.il  III  as  eiven  1“')  Holti  branches  were  lilted 
smiiill.ir-eoiislv 


2.2  Axial  Shock  Velocity 

Figure  9  shows  the  axial  velocity  of  the  shock  waves 
generated  by  nanosecond  pulses  at  two  energy  levels. 
The  positive  axis  corresponds  to  the  proximal  position, 
the  side  of  the  plasma  closest  to  the  incoming  laser  pulse. 
It  appears  that  the  shock  velocities  at  the  proximal  side 
are  higher  than  those  at  the  distal  side.  The  difference, 
however,  is  still  within  the  standard  deviation  of  the 
measurements.  It  is  possible  that  the  shock  velocity  and 
thus  the  pressure  is  higher  at  the  proximal  side.  This  could 
be  explained  by  the  way  the  plasma  expands  during  the 
laser  pulse;  the  proximal  side  practically  shielding  the 
distal  side  [28,  29].  In  the  theoretical  curves  presented 
here,  the  two  branches  have  been  fitted  separately. 

The  difficulty  in  the  analysis  of  the  axial  velocities  is 
that  the  center  of  the  plasma,  as  we  mentioned  earlier, 
cannot  be  determined  unambiguously.  Taking  a  diifcrent 
point  along  the  axial  direction  as  center  would  affect 
the  value  of  the  constant  (C)  considerably.  It  would  not, 
however,  affect  the  calculated  pressure.  This  is  because 
the  two-beam  method  gives  an  absolute  magnitude  of  the 
shock  velocity. 

Picosecond-generated  plasma  showed  no  discernible 
shock  waves  in  the  axial  direction.  Recent  experiments 
by  Vogel  and  co-workers  using  fast  photography  showed 
that  the  picosecond-generated  plasma  is  not  cotnpietely 
homogeneous  and  that,  at  the  proximal  side,  it  consists 
of  a  large  number  of  microplasmas  [30].  Figure  10  shows 
the  plasma  and  the  shock  wave  generated  by  a  3.6-mJ, 
30-ps  laser  pulse  30ns  after  optical  breakdown  [30].  In 
the  radial  direction  the  shock  waves  combine  to  create 
a  niea-.iir.ible  eHect.  while  in  the  axial  direction.  the\  aie 
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Fig.  10.  Picosecond-generated 
plasma  by  a  3.6-mJ  laser  pulse 
showing  a  large  number  of 
microplasma.  The  photograph 
was  taken  30  ns  after  the 
optical  breakdown.  Note  the 
appearance  of  the  shock  wave 
in  the  radial  direction  only. 

The  arrow  indicates  the  dir^on 
of  the  cxdting  laser  pulse.  The 
bar  equals  100  pnu  The 
photograph  was  taken  by 
Dr.  A.  Vogel  of  H.  Wacker 
Laboratory,  Munich  and  is 
reproduced  by  permission 


reflected  on  the  surface  of  the  neighboring  microplasmas 
and  dissipate.  The  large  impedance  difference  between 
the  liquid  and  gaseous  phases  accounts  for  the  efficient 
reflection  and  scattering  of  the  shock  waves,  so  that  no 
shock  waves  appear  in  the  axial  direction. 


3.  Discussion 

The  shock  wave  pressure  can  be  calculated  from  the 
shock  and  particle  velocities  and  the  jump  condition.  If 
(5)  and  (7)  are  substituted  into  (1),  we  derive  an  expres¬ 
sion  for  the  pressure  as  a  function  of  the  distance  from 
the  center  of  the  plasma 


Figures  11  and  12  show  the  calculated  shock  wave  pres¬ 
sures  in  the  radial  and  axial  directions,  respectively.  The 
present  work  seems  to  contradict  measurements  of  pres¬ 
sure  waves  using  transducers,  which  show  that  the  pres¬ 
sure  for  spherically  expanding  waves  is  proportional  to 
I  /r  [16]  rather  than  to  1  /r^.  It  should  be  noted,  however, 
that  the  pressure  calculated  from  our  velocity  measure¬ 
ments  relates  only  to  the  shock  waves.  The  shock  pressure 
decreases  rapidly  so  that  at  distances  over  500  pm  it  is 
practically  zero.  The  far  field  pressure  computed  from  the 
shock  velocity  is  not  accurate  because  the  velocity  of  the 
pressure  wave  is  very  close  tc>  the  velocity  of  .sound  and 
the  resolution  of  our  technic] i:  is  not  sufficient  to  allow 
for  an  accurate  measurement  c  .  the  velocity  in  this  range. 
Consequently,  the  pressure  in  the  far  field  is  underesti¬ 
mated.  Equation  (8)  is  valid  only  in  the  near  field  where 
the  pressure  waves  generated  by  the  optical  breakdown 
are  manifested  as  shock  waves. 

The  optical  measurements,  taken  together  with  the 
tran.sducer  measurements,  seem  to  support  the  concept 
that  the  plasma  initially  generates  a  shock  wave,  the 
pressure  of  which  is  proportional  to  l/r‘.  Ihe  shock 
wave  dissip'tcs  rapidlv  and  continuous  as  an  acoustic 
wave  that  decays  as  l/r.  T  he  dilTcrcncc  in  (lie  scaling 
l.i«s  t(ir  the  two  t\pc''  ol  waves  (snpcr.sonic  vs  sonicl 
rc'-iilis  Ironi  llic  anenualion  nu'chanisnis  In  ihc  sonic 
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RADIAL  DISTANCE  FROM  THE  CENTER  OF  THE  PLASMA  (mm) 

Fig.  11.  Shock  peak  pressure  vs.  radial  distance  calculated  from  (8). 
The  three  curv'es  correspond  to  the  three  sets  of  measurements  of 
the  shock  velocities  shown  in  Tips.  7  and  8 
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AXIAL  DISTANCE  FROM  THE  CENTER  OF  THE  PIJV.SMA  (mm) 

Fig-  12.  Shock  peak  pressure  vs.  axial  distance  calculated  from  (8) 
The  two  curves  correspond  to  the  two  sets  of  measurenienls  of  the 
shock  velocities  shown  in  l  ip.  9.  The  arrow  indicates  the  direction 
of  the  excilinp  laser  pulse 


region,  the  energy  of  the  wave  is  conserved  provided  that 
the  medium  does  not  absorb  the  sound.  For  a  spherically 
expanding  wave,  the  average  kinetic  energy  within  a  thin 
layer  at  ti  distance  r  fiom  the  center  ivf  the  plasma  can 
be  written  .is 
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tquation  (9)  leads  lo  a  pressure  that  is  proportional  to 
1/r.  A  shock  wave,  on  the  other  hand,  involves  irre¬ 
versible  processes,  wherein  the  kinetic  energy  within  the 
shock  front  is  converted  into  heat.  In  this  case,  the  mo¬ 
mentum  but  not  the  kinetic  energy  is  con  :erved.  Foi  a 
spherically  expanding  shock  wave  the  pressure  is  pro¬ 
portional  to  1/r^.  The  existence  of  two  different  math¬ 
ematical  expressions  for  the  pressure  imphes  that  the 
pressure  in  the  near  field  cannot  be  dedu«d  from  pres¬ 
sure  measurement,  in  the  far  field  unless,  of  course,  the 
proportionahty  constant  in  the  two  expressions  is  the 
same. 

It  can  also  be  instructive  to  compare  our  results  (8) 
with  two  theories  developed  for  underwater  explosions 
[31],  namely  Kirkwood-Bethe  and  Penney-Dasgupta.  Of 
the  two,  Penney-Dasgupta  theory  provides  a  simple  ex¬ 
pression  for  the  maximum  pressure  which  is  valid  in  the 
proximity  of  the  explosion, 

P  =  G(ao/r)  exp(2no/r)  (10) 

where  oq  is  the  linear  dimension  '‘f  the  explosive  and  G  a 
constant.  We  can  fit  (10)  to  our  measurements  taking  as  oq 
the  linear  dimension  of  the  plasma.  The  fit  is  good  except 
in  the  far  field  where  (10)  predicts  higher  pressures.  We 
know,  however,  that  our  determination  of  the  pressure  in 
the  far  field  is  not  accurate,  as  it  has  been  discus.sed  in 
the  experimental  section.  We  can  approximate  (10)  with 
the  first  two  terms  of  a  power  series  expansion,  since  ooA 
is  less  than  one  in  our  measurements. 

P  =G{ao/r)+2G(ao/rf  +  ...  .  (11) 

The  first  term  corresponds  to  the  acoustic  wave  that 
remains  after  the  dissipation  of  the  shock  wave.  The 
second  term  becomes  doniinant  in  the  near  field  and 
corresponds  to  the  shock  wave.  An  important  implication 
of  the  Penney-Dasgupta  expression  is  that  it  may  provide 
a  way  to  extrapolate  the  shock  pressure  in  the  near  field 
from  acoustic  pressure  measurements  in  the  far  field. 

Underwater  shock  wave  research  has  established  that 
the  shock  pressure  scales  as  the  linear  dimension  of  the 
explosive  [31].  The  fact  that  (10)  fits  our  measurements, 
when  oo  is  taken  as  the  linear  size  of  the  plasma,  indi¬ 
cates  that  a  similar  scaling  law  may  operate  in  the  case 
I'f  plasma-induced  shock  waves.  In  fact,  some  of  our 
measurements  in  the  radial  direction  are  consistent  with 
a  similar  scaling  law.  The  shock  pressures,  expressed  as 
a  function  of  normalized  distance  from  the  center  of  the 
plasma  r/ck).  arc  the  same.  The  difficulty  in  establishing 
this  scaling  law  unambiguously  in  our  experiments  is  re¬ 
lated  to  the  fluctuations  in  both  the  size  of  the  plasma 
and  the  shock  velocity.  In  our  measurements  the  shoek 
velocity  is  the  average  velocity  of  many  separate  shock 
waves  while  oo  is  the  smallest  linear  dimension.  What 
may  be  needed  is  the  simultaneous  measurement  of  the 
plasma  size  and  the  shock  velocity. 

4.  Conclusions 

\\.'c  luive  Used  a  iwn-licam  picihc  lo  meaMire  Ilie  shock 
vdociu  around  a  la^ci-m-iiiccd  pl.isni.i.  \k’e  have  demon 


strated  that  the  technique  is  accurate  enough  to  measure 
the  shock  velocity  and  hence  determine  the  shock  maxi¬ 
mum  pressure  non-invasively.  We  have  used  the  principle 
of  conservation  of  total  momentum  of  the  shock  front  to 
derive  expressions  for  the  shock  velocity,  particle  veloc¬ 
ity  and  pressure  vs.  the  distance  from  the  center  of  the 
plasma.  Our  work  supports  the  concept  that  the  expand¬ 
ing  plasma  initially  induces  a  shock  wave  which  dissipates 
rapidly  into  an  acoustic  w.ave.  In  addition,  we  have  found 
that  the  range  of  the  shock  wave  is  limited  compared  to 
the  range  of  acor:tic  waves.  Thus,  tissue  damage  caused 
by  shock  waves  should  also  have  a  limited  range. 
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Abstract  We  investigate  shock  waves  generated  by  excimer  laser  ablation  of  sheet  polyimide 
confined  in  water.  The  velocities  of  the  ablation-induced  pressure  waves  in  the  water  are 
determined  by  an  optical  probe  system.  We  measure  supersonic  velocities  up  to  a  few  hundred 
microns  away  from  the  irradiated  surface,  indicating  the  formation  of  shock  waves.  We  use 
these  velocities  to  calculate  the  corresponding  pressures.  They  are  already  in  the  kbar  range  at 
fluences  comparable  to  the  threshold  of  ablation.  The  shock  pressure  varies  as  the  square  root 
of  the  incident  laser  fluence,  a  behavior  that  is  explained  by  the  rapid  heating  of  the  confined 
gaseous  products  of  ablation. 

The  initially  planar  shock  waves  propagate,  become  spherical,  and  decay  within  a  few 
hundred  microns  in  the  surrounding  water  to  acoustic  waves.  During  spherical  expansion  the 
shock  pressure  drops  as  the  inverse  of  the  square  of  the  propagation  distance. 

The  shock  waves  generated  may  be  relevant  in  explaining  photoacoustic  damage  observed  in 
biological  tissue  after  excimer-ablation  at  corresponding  irradiances.  They  may  also  be  important 
in  material  processing  applications  of  excimer  laser  ablation  of  polymers  as  they  can  lead  to 
plastic  deformation. 

PACS:  81.60.Jw,  43.35.Cg,  87 


The  ability  of  the  ArF  excimer  laser  to  make  high  preci¬ 
sion  cuts  with  thermal  damage  extending  less  than  1  pm 
from  the  edge  of  the  incision  is  one  of  the  reasons  for 
the  interest  in  this  laser  for  corneal  and  keratorefractive 
surgery.  However  recent  studies  of  ArF  laser  ablation  of 
skin  in  air  have  shown  tissue  damage  up  to  200  pm  below 
the  bottom  of  the  incision  [1].  Similar  results  have  been 
found  earlier  in  cornea  [2,  3].  For  confined  ablation  of 
skin  in  water  the  damage  zone  was  significantly  larger, 
up  to  300  pm  [4].  Although  these  large  damage  zones 
have  all  been  attributed  to  photoacoustic  effects,  a  clear 
understanding  of  the  underlying  damage  mechanism  is 
missing.  In  particular  it  is  unclear  whether  a  shock  wave 
is  formed  or  if  the  damage  is  due  to  strong  acoustic 
waves. 

A  shock  wave  differs  qualitatively  from  an  acoustic 
pressure  transient.  It  is  characterized  by  a  supersonic 
propagation  speed  and  a  very  steep  front  that  for  all 
practical  purposes  can  be  treated  as  a  discontinuity  in 
pressure,  density,  particle  velocity  and  internal  energy. 
At  the  shock  front  the  pressure  rises  a  few  kbars  within 


a  distance  of  a  few  nanometers,  resulting  in  pressure 
transients  with  risetimes  of  a  few  picoseconds  [5].  A  shock 
also  changes  the  thermodynamic  state  of  the  medium  in 
which  it  propagates  in  a  different  way  than  an  adiabatic 
process  does. 

Acoustic  pressure  transients  of  several  hundred  bars 
with  risetimes  of  a  few  nanoseconds  have  been  obser\’ed 
in  excimer  ablation  of  tissue  and  organic  polymers  using 
polyvinylidenefluoride  (PVDF)  piezoelectric  transducers, 
[6,  7].  However  these  transducers,  because  of  their  limited 
bandwidth,  do  not  allow  detection  of  the  presence  or 
measurement  of  the  peak  pressure  of  shock  waves.  To 
date,  no  direct  detection  of  shock  waves  in  connection 
with  UV  excimer  ablation  has  been  reported. 

In  order  to  study  the  acoustic  effects  associated  with 
UV  laser  ablation  we  investigated  the  ablation  of  a  thin 
(25  pm)  sheet  of  polyimide  immersed  in  distilled  water. 
To  ablate  the  targets  we  use  a  XeCl  excimer  laser  emit¬ 
ting  at  308  nm,  as  the  high  optical  absorption  coeffiecient 
oo  of  polyimide  at  that  wavelength,  oo  =  10^cm“'  [8,  9], 
makes  it  a  reasonable  model  for  the  ablation  of  tissue 
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with  193  nm  radiation;  at  that  wavelength  the  absorp¬ 
tion  coefficient  of  the  cornea  has  been  measured  to  be 
2700cm-‘  [10], 

Using  an  optical  probe  system  we  measure  the  ve¬ 
locities  of  the  laser-induced  stress  waves  in  the  water  as 
a  function  of  incident  laser  fluence.  We  measure  super¬ 
sonic  velocities,  which  indicates  the  formation  of  shock 
waves,  at  fluences  close  to  the  threshold  of  ablation. 
While  it  would  be  desirable  to  analyze  the  generated 
shocks  directly  inside  the  polyimide,  the  optical  prop¬ 
erties  of  our  targets  do  not  allow  this  with  our  setup. 
However,  assuming  an  acoustic  reflection  at  the  interface 
between  the  two  materials,  the  characterization  of  the 
shock  waves  in  the  water  gives  physical  insight  into  the 
ablation  of  polyimide.  Our  work  also  suggests  that  the 
effects  of  laser-induced  stress  waves  may  be  important  in 
patterning  microelectronic  materials  by  UV-ablation. 

In  our  experiments  the  polyimide  sheet  is  completely 
immersed  in  water,  thus  increasing  photoacoustic  effects 
by  confinement  of  the  ablated  material.  Because  the  equa¬ 
tion  of  state  of  water  is  known,  using  jump  conditions  at 
the  shock  front  it  is  straightforward  to  calculate  the  peak 
pressures  corresponding  to  the  measured  shock  velocities. 
We  observe  the  generation  of  strong  (multi-kbar)  shocks 
in  water.  We  determine  how  these  pressures  decay  during 
the  expansion  of  the  shocks  in  the  water.  In  order  to  find 
the  thresholds  for  the  ablation  confined  in  water  and 
verify  that  the  logarithmic  dependence  of  etch-depth  on 
fluence  observed  for  ablation  in  air  applies  [9, 11],  we  also 
measure  the  etch-depth  per  pulse  of  the  immersed  poly¬ 
imide.  At  low  fluences  we  find  a  logarithmic  dependence 
on  the  incident  laser  fluence. 

Our  observations  are  explained  by  a  simple  ideal  gas 
model.  We  show  that  shocks  develop  from  the  high 
pressure  transient  that  is  generated  by  the  expansion 
of  rapidly  heated  confined  gas. 


1  Experimental  Methods 

Our  experimental  setup  is  illustrated  schematically  in 
Fig.  1.  To  ablate  our  targets  we  use  a  XeCl  excimer 
laser  (Lambda  Physik  model  EMG  103  MSC)  that  emits 
pulses  of  a  duration  of  20  ns  at  a  wavelength  of  308  nm 
[12].  A  circular  aperture  of  a  diameter  of  5  mm  selects 
a  uniform  part  of  the  laser  beam  that  passes  various 
attenuation  filters  and  an  aspheric  lens  with  a  focal  length 
of  18  mm.  The  beam  is  focused  through  a  quartz  window 
into  a  cuvette  that  is  filled  with  distilled  water  prior  to 
the  experiment.  During  the  experiment  the  water  is  not 
changed. 

As  targets  we  use  sheet  polyimide  of  a  thickness  of 
25  pm  immersed  in  water.  It  is  clear  that  the  confinement 
of  the  target  by  liquid  has  an  amplifying  influence  on  the 
acoustic  effects  unter  investigation.  To  keep  that  influ¬ 
ence  constant  the  irradiated  site  was  always  more  than 
10  mm  below  the  water  surface.  Indeed  it  has  been  shown 
in  experiments  on  laser  fragmentation  of  biliary  and  uri¬ 
nary  stones  performed  at  similar  irradiation  conditions 
as  discussed  here  that  the  amplitude  of  a  laser-generated 


Fig.  1.  Schematic  view  of  the  experimental  setup.  The  He-Ne  probe 
beams  are  separated  by  35  pm  and  can  be  translated  in  space 
without  changing  their  relative  spatial  separation.  The  target  T  is 
mounted  on  a  separate  translation  stage  and  can  be  positioned 
relative  to  the  probe  beams  with  an  accuracy  of  10  pm.  The  z- 
axis  is  the  axis  of  the  excimer  beam  EB  with  z  =  0  at  the  rear 
surface  of  the  target.  To  determine  penetration  of  the  target,  the 
transmitted  UV-radiation  is  collected  with  lens  L3  and  imaged 
onto  photodiode  D3.  C  cuvette  filled  with  water,  BS  beam  splitter, 
M  mirror,  L,  lenses,  £>/  detectors,  F|  neutral  density  filter,  Fz 
bandpass  filter,  A  aperture,  A/D  oigitizer 


stress  wave  does  not  depend  on  the  immersion  depth, 
provided  that  that  is  larger  than  3  mm  [13]. 

To  determine  the  velocity  of  shock  waves  propagating 
in  water  we  use  an  optical  probe  technique  that  was 
described  in  detail  elsewhere  [14].  We  focus  two  TEMoo 
Helium-Neon  probe  beams  into  our  cuvette  in  such  a  way 
that  they  are  parallel  to  each  other  and  to  the  surface  of 
the  target.  A  fast  photodiode  (EG&G  model  FND  100, 
risetime  <  1  ns)  detects  them  after  they  have  propagated 
through  an  imaging  lens  and  an  aperture. 

The  probe  beams  have  a  waist  of  cuo  =  7  iim  and  are 
laterally  separated  by  a  distance  of  about  35  pm.  Their 
axes  intersect  the  axis  of  the  ablating  beam,  which  is  per¬ 
pendicular  to  the  surface  of  the  target  and  denoted  as  the 
z-axis  in  Fig.  1.  The  intersection  between  the  back  surface 
of  the  target  and  the  z-axis  is  taken  as  the  origin  z  =  0. 
X-Z -translation  stages  allow  the  simultaneous  displace¬ 
ment  of  the  two  probe  beams  parallel  and  perpendicular 
to  the  surface  of  the  target.  Unless  otherwise  specified 
the  probe  beams  are  at  the  back  of  the  irradiated  targets. 

The  pressure  waves  propagating  in  the  water  create 
density  variations  and  therefore  transient  changes  in  the 
index  of  refraction.  As  they  traverse  the  probe  beams, 
these  beams  are  deflected  out  of  the  acceptance  angle  of 
the  detection  optics.  Because  the  changes  in  the  index 
of  refraction  of  the  water  are  small,  on  the  order  of  a 
few  percent  [15],  deflection  is  efficient  only  near  graz¬ 
ing  incidence,  so  that  once  the  pressure  wave  has  moved 
through  the  probe  beam  at  its  waist,  the  probe  beam  is 
again  imaged  onto  the  photodiode.  The  net  result  is  a 
decrease  of  the  photodiode  signal  (D2)  for  a  duration 
roughly  equal  to  the  transit  time  of  the  pressure  wave 
through  the  probe  beam,  about  ten  nanoseconds  in  our 
setup.  Because  both  probe  beams  produce  the  same  effect, 
the  photodiode  signal  consists  of  two  such  variations  in 
sequence.  The  signal  is  digitized  by  a  Tektronix  digitizer 
(model  AD 7912,  750 MHz  bandwidth)  and  sent  on  a 
GPIB  bus  to  an  IBM-AT-clone  computer  for  analysis.  A 
computer  program  locates  the  two  peaks  of  the  photodi- 
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ode  signal  and  allows  us  to  determine  their  separation  to 
an  accuracy  of  0.5  ns.  That  time  is  the  transit  time  of  the 
pressure  wave  through  the  axes  of  our  two  probe  beams. 
At  sufficiently  large  distances  from  the  target,  the  pres¬ 
sure  waves  become  sonic  and  propagate  at  the  velocity 
of  sound  in  water  of  20°  C,  Cs  =  1.483  km/s.  Measuring 
the  transit  time  far  away  from  the  target  thus  allows  us 
to  determine  the  distance  between  the  two  probe  beams. 

In  short,  our  setup  allows  the  non-invasive  determina¬ 
tion  of  shock  velocities  in  water  to  an  accuracy  of  5%. 
The  spatial  resolution  is  equal  to  the  axial  separation  of 
the  two  probe  beams,  ~  35  pm.  In  addition,  we  can  mon¬ 
itor  the  dynamics  of  gas  bubbles  created  during  ablation 
by  locating  the  probe  beams  in  front  of  the  target  where 
the  ablation  takes  place.  As  cavitation  bubbles  typically 
expand  at  velocities  of  20-100  m/s,  their  signals  can  eas¬ 
ily  be  distinguished  from  those  created  by  the  more  than 
fifteen  times  faster  pressure  waves. 

When  measuring  the  shock  velocities,  we  changed  the 
targets  frequently  in  order  to  eliminate  the  effects  of 
thinning  of  the  target  by  ablation;  at  most  three  laser 
pulses  were  fired  onto  the  same  surface  spot.  Every  time 
a  fresh  spot  was  irradiated,  the  target  was  repositioned 
carefully  relative  to  the  probe  beams.  This  is  imperative, 
as  the  shock  waves  transmitted  into  the  water  attenuate 
significantly  within  a  few  tens  of  microns  because  of 
geometrical  expansion.  Ideally  the  axis  of  the  ablating 
beam  and  the  two  waists  of  the  probe  beams  should  lie 
in  the  same  plane.  The  axis  of  the  ablating  beam  should 
also  intersect  the  surface  of  the  target  and  the  axes  of 
the  two  probe  beams  at  a  right  angle.  Moving  the  target 
through  the  probe  beams  and  simultaneously  measuring 
the  transmitted  light  on  the  photodiode  (D2)  allows  us 
to  position  the  target  to  an  accuracy  of  10  pm. 

To  determine  the  etch  rate  of  polyimide  immersed  in 
water  we  measured  the  number  of  pulses  to  penetrate. 
We  added  a  lens  (L3)  that  images  the  UV-light  transmit¬ 
ted  through  the  sample  onto  another  photodiode  (D3) 
placed  behind  an  UV-bandpass  filter  to  the  above  setup. 
Since  the  optical  absorption  coefficient  oq  of  polyimide  at 
308  nm  is  very  high,  there  is  no  signal  from  that  diode  un¬ 
less  the  target  has  been  nearly  penetrated.  It  is  therefore 
relatively  easy  to  count  the  number  of  pulses  necessary 
to  drill  through  a  target  sheet  of  polyimide.  Dividing 
its  thickness  by  that  number  gives  the  etch-depth  per 
pulse.  We  took  data  from  fluences  of  ~  lOOJ/cm^  down 
to  ~  150mJ/cm^.  Once  the  number  of  pulses  to  drill 
through  became  larger  than  5000,  corresponding  to  an 
etch-depth  per  pulse  of  less  than  5  nm,  we  assumed  that 
the  fluence  was  below  the  ablation  threshold  and  stopped 
the  data  acquisition. 

The  fluence  F  of  the  ablating  pulse  at  the  surface  of 
the  target  is  calculated  from  measurements  of  the  pulse 
energy  and  the  spot  size.  We  reduce  the  maximum  pulse 
energy  of  about  4mJ  by  a  set  of  attenuating  filters  (FI) 
and  determine  the  energy  of  the  ablating  pulse  by  mea¬ 
suring  a  split-off  part  of  the  beam  with  a  pyroelectric 
detector  (Dl)  (Molectron  model  J3-09).  To  obtain  the 
energy  deposited  into  the  focal  spot,  we  multiply  the 
measured  energy  with  a  transmission  factor  that  corrects 
for  the  attenuation  of  the  beam  because  of  reflection  at 


quartz-air  and  quartz-waicr  interfaces.  The  area  of  the 
irradiated  spot  is  determined  by  microscopic  examina¬ 
tion.  It  is  elliptical  with  principal  axes  of  90  and  62  pm 
and  an  area  A  =  4.4  x  10^^  cm*.  We  calculate  the  fluence 
F  =  Etn/A  by  assuming  uniform  illumination  of  the  fo¬ 
cal  spot  with  a  laser  pulse  of  energy  We  work  with 
fluences  up  to  90J/cm‘  corresponding  to  irradiances  up 
to  4.5GW/cm^. 

2  Results 

2.]  Shock  Wave  Measurements 

Below  a  fluence  of  8mJ/cm^  we  do  not  obtain  any  sig¬ 
nal  with  our  optical  probe-system,  neither  in  front  nor 
in  back  of  the  irradiated  targets.  Starting  at  that  fluence 
value  small  gas  bubbles  form  at  the  irradiated  surface. 
As  they  scatter  a  probe  beam  that  illuminates  the  front 
surface  at  grazing  incidence,  we  can  detect  them  with 
our  setup.  They  may  consist  of  the  gases  that  have  been 
observed  to  emanate  from  below  the  surface  of  polyimide 
during  ablation  in  air  at  very  low  fluences  [16]  or  simply 
of  water  vapor  resulting  from  conductive  heating  of  the 
water  by  the  irradiated  surface  of  the  target.  At  these 
low  fluences  the  generated  microcavities  have  lifetimes  of 
about  10  ps.  The  acoustic  transients  that  they  generate 
are  below  our  detection  limit. 

At  a  fluence  of  Fo  =  50mJ/cm^  we  begin  detecting 
pressure  waves  in  front  of  the  target.  They  propagate  at  a 
speed  about  5%  or  75  m/s  above  the  speed  of  sound.  The 
corresponding  pressure  transients  that  are  launched  into 
the  polyimide  are  still  too  small  to  be  detected  by  our 
probe  system  at  the  back  side  of  the  target.  Because  Fq 
is  about  equal  to  the  threshold  of  ablation  determined  in 
air  by  photoacoustic  techniques  (45mJ/cm^  [11,  17]),  we 
conclude  that  Fo  is  the  threshold  for  confined  ablation, 
even  though  it  is  smaller  than  the  value  fs  determined 
from  the  extrapolation  of  the  etch-depth  data  below  (cf. 
Sect.  2.2). 

For  fluences  between  130mJ/cm^  and  about  KX)J/cm^ 
we  monitor  the  strength  of  the  transmitted  shock  waves. 
During  spherical  expansion  in  the  water,  their  velocities 
drop  to  the  speed  of  sound  after  propagation  over  a 
distance  of  a  few  hundred  microns. 

We  calculate  the  peak  pressures  of  the  shock  waves 
propagating  in  water  from  measured  shock  velocities  by 
using  the  jump  condition, 

f’s  =  eoF'sUp,  (1) 

where  Fs  is  the  shock  pressure,  qo  =  0.987  g/cm^  the 
density  of  water  at  20°  C  and  and  Wp  the  shock  and 
particle  velocities,  t/s  and  Up  are  further  related  by  the 
equation  of  state.  To  first  approximation 

Us=^a  +  hup,  (2) 

with  the  values  a  =  1.483  km/s  and  h  =  2.07  for  water 
[15]. 

Shock  pressures  determined  from  the  measured  shock 
velocities  are  presented  in  Fig.  2a,  b.  Figure  2a  shows 
the  dependence  of  the  shock  pre.ssure  PJzo.F)  on  the 
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Fig.  2a,  b.  Shock  pressures  of  the  waves  transmitted  into  the  water 
at  the  rear  surface  of  25  pm-thick  samples  were  calculated  from 
measured  velocities  using  (1-2).  a  Shock  pressure  P,(2o,f)  vs  flu- 
enoe  F  at  an  =  25  pm  (o),  50  pm  away  from  the  irradiated  surface. 
The  dashed  line  is  a  least  square  ht  of  the  data  in  double-logarithmic 
representation.  It  has  a  slope  of  0.52.  The  solid  lines  are  predictions 
of  the  thermal  model  explained  in  the  text.  They  were  calculated 
with  (6)  assuming  the  expanding  gas  to  be  ideal  with  a  molar  spe¬ 
cific  heat  at  constant  volume,  Cy  =  {f  /2)R,  where  /  is  either  5  or 
12.  b  Shock  pressure  P,(z,f,)  vs  distance  z  from  the  rear  surface 
of  the  target  for  two  fixed  laser  fluences  F\  =  16.5  J/cm^  (•)  and 
Fi  =  0.7  J/cm^  (o).  The  solid  lines  are  fitted  curves  that  assume  a 
1/(2  —  r,)^  decay  according  to  (3).  Such  a  decay  is  expected  from 
momentum  conservation  of  the  shock  front 


fluence  F  of  the  ablating  pulse  at  the  fixed  distance 
2o  =  25  pm  away  from  the  rear  surface  of  the  target. 
The  pressure  data  plotted  in  a  double-logarithmic  plot 
are  well  fitted  by  a  straight  line  (dashed),  implying  a 
power  law  dependence  between  the  shock  pressure  and 
the  laser  fluence.  The  least  square  fit  gives  a  slope  of  0.52. 
At  250mJ/cm^,  5  times  the  threshold  of  ablation,  which 
is  in  the  typical  working  range  of  practical  applications, 
the  fit  indicates  a  pressure  Ps  of  905  bar.  The  two  solid 
lines  in  Fig.  2a  are  theoretical  predictions  from  the  ideal 
gas  model  discussed  in  Sect.  3. 1 .  They  were  calculated  by 
assuming  the  gas  to  have  a  molar  specific  heat  at  constant 
volume  Cy  of  either  5/2  R  or  6  R,  R  being  the  universal 
gas  constant. 

In  Fig.  2b  we  show  the  decay  of  the  shock  pressure 
with  distance  z  from  the  rear  surface  of  the  target  at 
z  =  0  as  the  shock  wave  propagates  into  the  surround¬ 
ing  water  for  the  two  fluences  F[  =  16.5  J/cm^  and 
F2  =  0.7J/cm^.  At  large  enough  distances  z  the  shock 
waves  expand  spherically,  and,  because  of  momentum 
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conservation,  their  pressure  /\  varies  inversely  with  the 
square  of  the  distance  from  the  center  of  the  shock  [14]. 
Accordingly  we  fit  the  data  by ; 

Ps(z,F,)  =  F,(^^^^  ,  (3) 

where  Pj  and  2,,  i  =  1,2,  are  the  two  parameters  of  that  fit. 
We  obtain  numerically  Pi  =  8.3kbar  and  P2  =  1.9kbar. 
These  values  agree  well  with  the  data  from  Fig.  2a,  which 
is  expected,  as  P,  =  Ps(zo,P,).  The  other  parameter  of 
(3),  z„  indicates  where  on  the  z-axis  the  center  of  the 
spherically  expanding  shock  lies.  If  the  ablation-induced 
shock  were  always  spherical,  we  would  expect  the  center 
of  shock  to  lie  on  the  irradiated  front  surface  of  the  target, 
Zj  X  —25  pm.  If  it  were  planar,  we  would  have  z,  =  —  00. 
The  fitted  values,  zj  =  —  88  pm  and  zj  =  —  75  pm,  are  in 
between,  indicating  that  the  generated  shocks  are  initially 
planar  and  become  spherical  only  after  propagating  a  few 
tens  of  microns. 


2.2  Etch-Rate  Measurements 

At  fluences  equal  to  several  times  the  ablation  threshold, 
photo-etching  of  polyimide  at  the  XeCl  laser  wavelength 
of  308  nm  is  generally  well  described  by  a  blow-off  model. 
In  that  model  one  assumes  that  the  radiation  attenuates 
exponentially  with  distance  from  the  surface  into  the 
target  according  to  some  effective  absorption  coefficient 
Oe,  and  that  all  material  that  absorbs  an  energy  density 
>  oteFs,  where  Fj  is  a  threshold  fluence  for  significant 
photo-etching,  is  ablated.  The  resulting  etch-depth  d  is 
given  by  [9,  18,  19]: 

d  =  a;'ln(^~^.  (4) 

Our  etch-data,  shown  in  Fig.  3  in  a  semi-logarithmic  rep¬ 
resentation,  follows  (4)  for  fluences  below  about  20J/cm^ 
with  the  fitted  parameter  Fs  =  150mJ/cm^  and  0^  = 
3.6  X  lO^cm"’.  These  values  are  similar  to  those  obtained 
by  other  workers  for  photo-etching  of  polyimide  in  air 
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Fig.  3.  Etch-depth  d  vs  fluence  F  for  confined  ablation  of  poly¬ 
imide.  The  solid  line  is  a  fit  to  the  data  by  t4)  for  fluences  less 
than  20J/cm^.  The  sharp  increase  of  etch-depth  at  that  fluence  is 
accompanied  by  permanent  mechanical  deformation  of  the  target 
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Table  1,  Kflective  absorption  coeflicient  ct^  and  threshold  of  ablation 
fs,  derived  from  published  data  by  assuming  a  blowolT  model 
according  to  (4).  When  the  data  deviated  from  that  dependence  at 
high  fiuences,  only  the  low-flucnce  part  was  used  for  the  fit.  Note 
that  dc  is  smaller  than  the  low  fluence  optical  absorption  coefficient 
do  =  10* cm"',  and  that  is  larger  than  the  threshold  fiuences 
determined  photoacoustically,  f'o  =  45mJ/cm’ 
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(cf.  Table  1),  indicating  that,  at  least  for  low  fiuences,  the 
etch-depth  of  polyimide  is  not  influenced  by  the  confine¬ 
ment  of  the  target. 

At  fiuences  above  20J/cm^  our  etch-depth  per  pulse 
becomes  larger  than  what  (4)  predicts  and  the  micro¬ 
scopic  appearance  of  the  ablated  region  is  different  from 
that  resulting  from  irradiation  a  lower  fiuences:  At  the 
front  surface  the  drilled  holes  are  now  surrounded  by 
a  funnel  shaped  region  of  modified  surface  texture  that 
indents  slightly  into  the  target.  Width  and  depth  of  that 
region  increases  with  fluence  and  amount  to  about  30  pm 
and  9  pm  at  90J/cm^;  we  attribute  this  to  the  perma¬ 
nent  mechanical  deformation  of  the  target  by  the  strong 
pressure  transients. 


3  Discussion 

3.1  Shock-Pressure 


To  calculate  the  generated  shock  pressures  we  adapt  an 
ideal  gas  model  used  by  other  workers  to  explain  shock 
generation  from  confined  ablation  of  metals  [20,  21].  In 
that  picture  it  is  assumed  that  a  constant  mass  of  an  ideal 
gas  enclosed  by  two  pistons  is  heated  by  a  laser  pulse.  If 
the  expansion  of  the  gas  is  violent  enough  to  drive  shock 
waves  into  the  pistons,  the  two  confining  surfaces  of  the 
expanding  gas  move  at  the  particle  velocities  of  the  piston 
material  through  which  these  shock  waves  have  passed. 
The  correct  boundary  conditions  of  the  expansion  of  the 
gas  thus  follows  from  the  relationship  between  shock 
pressure  and  particle  velocity  of  the  confining  materials. 
The  resulting  peak  pressure  Pp  is  given  in  [20]  as: 


Pp(f )  =  0.95  ■ 


ZF 

(I  +//2)t 


1/2 


(5) 


where  P  is  the  laser  fluence  absorbed  by  the  gas,  r 
the  duration  of  the  laser  pulse,  Z  an  effective  shock 
impedance  and  /  =  2c\/R,  Cy  being  the  molar  specific 
heat  of  the  ideal  gas  and  P  the  universal  gas  constant. 


The  factor  0.95  in  (5)  was  derived  in  a  dimensionle.ss 
numerical  analysis  of  the  gas  expansion  that  assumed 
heating  by  absorption  of  a  laser  pulse  of  Gaussion  time 
profile. 

Z  =  Z|Z2/(Z|  +  Z2)  is  an  efi'cctive  impedance  that 
combines  the  shock  impedances  Z,  of  the  two  materials, 
i.e.  the  respective  products  of  density  and  shock  velocity. 
We  shall  approximate  them  below  by  Z|  =  30.7  bar/(m/s) 
for  polyimide  [22]  and  Zj  =  14.6  bar/(m/s)  for  water,  the 
products  of  density  and  adiabatic  sound  speed.  This  is 
common  in  the  approximate  theory  of  traveling  waves  of 
finite  amplitudes  [23]. 

Note  that  (5)  implies  that  a  constant  fraction  of  the 
energy  that  heats  the  gas  is  converted  into  acoustic  en¬ 
ergy.  That  fraction  depends  only  on  the  impedances  in¬ 
volved  and,  through  /,  on  the  specific  heat  of  the  gas. 
This  follows  because  the  intensity  of  an  acoustic  wave  is 
proportional  to  the  square  of  its  pressure. 

In  thermal  coupling  measurements  it  has  been  found 
that  when  irradiating  polyimide  below  a  threshold  flu¬ 
ence  Fq,  all  the  absorbed  laser  energy  Ei„  is  converted 
into  heating  of  the  target  [17].  At  fiuences  above  that 
threshold  the  heat  deposited  in  the  target  remains  ap¬ 
proximately  constant,  and  the  excess  energy  —  AFo,  A 
being  the  area  of  the  irradiated  spot  and  Fq  the  thresh¬ 
old  fluence  of  ablation,  is  carried  away  by  the  ablated 
material  [17].  Consequently  we  modify  the  above  (5)  by 
replacing  f  by  F  —  fo  as  the  heat  source  of  the  confined 
gas. 

Because  of  the  mismatch  of  the  acoustic  impedances  at 
the  interface  between  polyimide  and  water,  the  pressure 
created  during  ablation  at  the  polyimide  surface  is  higher 
than  the  peak  pressure  of  the  shock  propagating  in  the 
water.  As  the  shock  passes  the  material  interface,  it  is 
only  partially  transmitted.  The  shock  pressure  is  reduced 
by  a  factor  Tac  =  2/(1  -I-Z1/Z2),  which  equals  0.65  in  our 
situation.  The  shock  pressure  Fs  of  the  waves  propagating 
into  the  water  at  the  rear  side  of  the  target  is  then  given 
by  the  ablation  pressure  Fp  from  (5)  multiplied  by  Tac- 
According  we  have: 

W)-0.95  .  (6) 

In  this  equation  the  number  /,  which  is  related  to  the 
molar  specific  heat  Cy  of  the  gas  by  /  =  2cy/R,  is  the 
only  parameter  that  is  not  well  known.  If  the  molecular 
vibrations  are  frozen,  which  is  a  good  approximation  at 
room  temperature,  /  equals  the  number  of  translational 
and  rotational  degrees  of  freedom  of  a  gas  molecule.  / 
becomes  larger  at  high  temperatures  where  a  considerable 
part  of  the  gas  molecules  are  vibrationally  excited.  The 
contribution  of  each  normal  vibration  of  each  molecule  to 
the  specific  heat  Cy  increases  with  temperature  T,  reaching 
the  value  k^T  asymptotically,  k^  being  the  Boltzmann 
constant.  Therefore  /  is  equal  to  or  larger  than  the 
number  of  tranlational  and  rotational  degrees  of  freedom 
of  a  molecule  and  smaller  than  that  number  increased  by 
twice  the  number  of  vibrational  degrees  of  freedom  [24]. 

The  principal  gaseous  ablation  products  of  polyimide 
in  air  are  CO,  CO2  and  H2O  and  elemental  carbon  [9,  25]. 
In  vacuum,  in  addition  IICN  and  ben/enc,  which  how- 
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ever  may  not  be  a  primary  ablation  product  [26],  have 
been  detected  [9,  25],  In  order  to  describe  the  gaseous 
ablation  products  by  a  simple  ideal  gas,  admittedly  a 
crude  approximation,  we  consider  below  the  two  limiting 
cases  /  =  5  and  /  =  12.  /  =  5  corresponds  to  a  diatomic 
or  linear  polyatomic  molecule  with  frozen  vibrations  and 
/  =  12  follows  for  a  gas  consisting  of  nonlinear  tri- 
atomic  molecules,  e.g.  H2O,  where  all  the  three  normal 
vibrations  are  contributing  the  maximum  value  R  to  the 
specific  heat.  As  the  principal  gaseous  ablation  products 
are  di-  and  tri-atomic  molecules  these  two  values  for  / 
give  the  extremal  predictions  of  the  model  for  the  pres¬ 
sure.  As  can  be  seen  from  (6)  they  differ  by  a  factor  of 
^/2. 

We  take  for  Fq,  the  threshold  of  ablation,  the  lowest 
fluence  where  our  probe  system  detects  the  emission  of 
pressure  waves,  Fq  =  50mJ/cm^  (cf.  Sect.  2.2),  and  we 
assume  a  pulse  duration  of  t  =  20  ns.  Note  that  the 
prediction  of  (6)  is  of  course  not  sensitive  to  the  exact 
value  of  Fq  at  high  ffuences. 

The  resulting  pressure  versus  fluence  dependences  cor¬ 
responding  to  /  =  5  and  /  =  12,  calculated  by  (6),  are 
plotted  in  Fig.  2a  as  solid  lines.  Given  the  simplicity  of  the 
model,  the  congruence  of  its  prediction  with  the  exper¬ 
imental  data  is  remarkable.  However  in  the  calculation 
several  simplifications  have  been  made  and  must  be  kept 
in  mind.  First  of  all,  the  calculation  assumes  a  planar 
shock.  The  measurement  shown  in  Fig.  2a  was  at  50  pm 
away  from  the  irradiated  surface.  That  distance  is  only 
a  little  bit  less  than  the  small  principal  diameter  of  our 
focal  spot,  62  pm,  and  edge  effects  have  probably  already 
started  to  weaken  the  shock.  Also  the  illumination  of 
the  focal  spot  is  not  perfectly  uniform.  Second,  the  mass 
of  the  removed  material  increases  with  fluence.  The  cor¬ 
responding  increase  of  the  separation  of  the  pistons  that 
confine  the  expanding  gas  in  the  model  description  is 
equal  to  the  etch  depth.  In  our  experiment  this  varies  be¬ 
tween  9%  and  26%  of  the  piston  separation  at  the  end  of 
the  laser  pulse,  which  is  predicted  by  the  model  as  PpX  jZ 
[21].  Increasing  the  piston  separation,  while  keeping  the 
total  internal  energy  constant,  reduces  the  pressure  in  the 
gas.  As  this  effect  is  neglected  in  the  model,  the  predicted 
pressures  should  slightly  overestimate  the  experimental 
data.  Finally,  as  mentioned  above,  we  assumed  in  the  cal¬ 
culation  an  ideal  gas  with  either  no  or  three  fully  excited 
normal  vibrations  per  molecule. 

The  pressure  of  the  ablation-induced  shock  waves 
during  spherical  expansion  in  the  water  varies  as  the 
reciprocal  of  the  square  of  the  distance  from  the  center 
of  the  shock.  This  has  been  observed  earlier  for  shock 
waves  induced  by  optical  breakdown  and  explained  by 
conservation  of  momentum  of  the  shock  front  [14].  Note 
however  that  before  spherical  expansion  becomes  impor¬ 
tant,  the  shock  propagates  as  a  plane  wave  for  a  distance 
that  increases  with  the  spot  size  of  the  ablating  beam. 

At  an  irradiance  of  ~  4.5  GW/cm^  we  generate  shocks 
of  ~  30kbar  inside  the  polyimide  using  pulses  of  4mJ. 
These  pressures  are  comparable  to  those  of  shocks  gen¬ 
erated  in  confined  laser  ablation  of  metal  foils.  There 
pressures  of  about  40  kbar  were  reported  for  irradiances 
of  a  few  GW/cm^  obtained  using  laser  pulses  of  a  dura¬ 


tion  of  30ns  and  pulse  energies  of  up  to  100  J  [21].  Note 
however  that  these  high-energy  laser  systems  illuminated 
focal  spots  of  6  mm  diameter,  producing  extended  planar 
shocks. 

3.2  Etch-Depth 

It  is  important  to  note  that  (4)  fails  to  describe  the  etch- 
depth  data  near  the  threshold  of  ablation  [16].  Using 
quartz  crystal  microbalances  it  has  been  shown  that  near 
the  threshold  of  ablation  the  amount  of  removed  material 
decreases  less  rapidly  with  fluence  than  predicted  by  (4) 
[27, 28].  It  was  concluded  that  photo-etching  of  polyimide 
at  308  nm  is  a  thermally  controlled  process,  with  ablation 
starting  at  surface  temperatures  of  about  850°  C  [28]. 

Therefore  (4),  although  a  good  approximation  at 
higher  ffuences,  becomes  an  inappropriate  description 
of  the  ablation  process  very  close  to  threshold.  It  fol¬ 
lows  that  Fs  and  0^,  the  two  parameters  of  the  blow-off 
model  that  are  calculated  by  fitting  etch-depth  to  (4), 
are  generally  different  from  the  actual  ablation  threshold 
Fq  and  the  optical  absorption  coefficient  oq.  This  is  con¬ 
firmed  by  the  values  in  Table  1 .  There,  for  all  experiments 
performed  with  nanosecond  pulses,  fs  is  larger  than  the 
threshold  of  ablation  found  photoacoustically  and  the 
effective  absorption  coefficient  is  smaller  than  oo,  the 
optical  absorption  coefficient  determined  by  transmission 
measurements  at  sub-ablative  ffuences.  Note  that  the  ex¬ 
perimental  data  obtained  by  the  different  workers  are 
consistent  in  that  the  product  FsO*  is  always  comparable 
to  and  larger  than  Fq(xq  =  4.5  kJ/cm^. 

At  ffuences  higher  than  20J/cm^  our  etch-depth  per 
pulse  data  deviate  from  (4).  We  observe  permanent  me¬ 
chanical  deformation  of  the  target  sheets,  and  the  micro¬ 
scopic  appearance  of  the  ablated  regions  changes.  Re¬ 
move  Sect. 3.3  from  the  manuscript:  Move  the  first  two 
sentence  of  Sect.  3.3  (“At  ffuences  higher  than  ...  the  ab¬ 
lated  regions  changes”)  as  a  separate  paragraph  to  the 
end  of  Sect.  3.2.  Delete  the  remaining  text  of  Sect.  3.3. 

4  Conclusions 

The  shock  waves  generated  by  ablation  of  polyimide 
confined  in  water  using  308  nm  excimer  laser  radiation 
have  been  studied.  The  measured  shock  pressures  are 
proportional  to  the  square  root  of  the  fluence  of  the 
ablating  pulse.  This  can  be  explained  by  a  simple  thermal 
model  that  describes  the  expansion  of  a  confined  ideal  gas 
subjected  to  rapid  heating.  According  to  the  model,  the 
ablation  pressure  depends  on  irradiance,  on  the  specific 
heat  of  the  gas  mixture  created  during  the  ablation  and  on 
the  acoustic  impedances  of  the  target  and  the  confining 
material. 

At  fluence  F  >  3(X)  mJ /cm^,  about  twice  the  threshold 
for  significant  photo-etching,  shock  pressures  of  more 
that  1  kbar  are  launched  into  the  confining  water.  Tissue 
ablation  by  excimer  lasers  is  usually  carried  out  at  similar 
irradiances.  It  is  therefore  likely  that  true  shocks  and 
not  just  strong  acoustic  waves  are  generated  in  tissue 
during  ablation  by  laser  pulses  of  a  duration  of  a  few 
nanoseconds  or  less. 
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A  steady-state  model  for  ablation  of  materials  by  strongly  absorbed  long-pulse  lasers  is 
developed.  Material  is  removed  by  a  combination  of  evaporation  and  ejection  of  liquid  due  to 
evaporation-induced  pressure  gradients  at  the  bottom  of  the  ablation  crater.  The  relative 
strength  of  the  two  processes  depends  on  the  details  of  the  liquid  flow.  An  explicit  expression 
for  the  drilling  velocity  is  derived.  The  model  explains  experimental  drilling-depth  and 
recoil  data.  It  is  applicable  to  materials  that  melt  and  vaporize. 


I.  INTRODUCTION 

Drilling  depth  and  target  recoil  provide  important  and 
easily  measurable  information  on  laser  ablation.  They  vary 
with  distance  between  the  surface  of  the  target  and  the 
laser  beam  waist.  For  pulses  of  several  100-/rs  duration,  the 
observed  variation  is  very  similar  for  tissue,  gelatin,  and 
metals,  suggesting  that  the  same  physical  mechanisms  are 
involved  when  ablating  these  materials.''^ 

Modeling  ablation  by  evaporation  jmd  ejection  is  not  a 
new  idea,  but  has  long  been  used  to  describe  laser  ablation 
of  metals^*  and,  more  recently,  of  biological  targets.*’’ 
Here,  I  will  show  how  a  model  based  on  such  an  approach 
indeed  explains  the  experimentally  observed  dependence  of 
drilling  depth  and  recoil  momentum  on  waist-surface  dis¬ 
tance.  The  concept  is  equally  well  applicable  to  metals  and 
biological  targets,  since  these  two  types  of  target  materials 
differ  primarily  in  the  physical  mechanism  that  determines 
the  depth  of  energy  deposition.  In  metals  it  is  heat  conduc¬ 
tion,  while  in  tissue  it  is  the  absorption  of  radiation.  The 
ablation  of  gelatin  by  pulsed  Er’''':YAG  laser  radiation, 
whose  wavelength,  2.94  fim,  corresponds  to  the  peak  of  the 
water  absorption  in  the  infrared,  will  be  used  to  illustrate 
the  model. 

The  steady-state  model  describes  ablation  of  materials 
by  strongly  absorbed  laser  radiation.  It  neglects  heat  con¬ 
duction  and  assumes  that  the  optical  penetration  depth  of 
the  radiation  is  small  compared  to  the  spot  diameter  of  the 
drilling  beam.  Attention  will  be  given  to  the  influence  of 
the  spot  size  on  the  ablation  mechanism.  The  drilling 
pulses  are  assumed  to  last  100  fis  or  longer,  corresponding 
to  pulsed  emission  of  free-running  solid-state  lasers  and 
axial-flow  CXI)2  lasers.  The  model  assumes  irradiances  well 
above  threshold,  so  that  a  strong  vaporization  develops 
instantaneously  at  the  target  surface.  The  radial  gradient  of 
the  evaporation-induced  pressure  on  the  target  surface 
causes  the  displacement  and  eventual  ejection  of  liquified 
material.  This  increases  the  drilling  eflSciency  and  affects 
the  recoil  of  the  target  induced  by  the  gaseous  and  liquid 
ejecta. 

Note  that  the  model  assumes  the  formation  of  a  liquid 
layer  at  the  surface  of  the  target  during  ablation.  It  there¬ 
fore  describes  the  ablation  of  materials  that  form  a  molten 
phase  before  vaporization.  It  does  not  apply  to  ablation  by 
photochemical  mechanisms  in  which  the  products  are 


formed  directly  by  decomposition  of  the  target,  for  exam- 
pie,  the  ablation  of  some  polymers  by  UV  light. 

First,  a  description  of  the  laser-induced  evaporation 
will  be  given,  assuming  the  formation  of  a  Knudsen  layer 
at  the  target  surface  and  taking  into  account  the  effect  of 
atmospheric  pressure.*  Then,  the  effect  of  radial  pressure 
gradients  on  a  thin  surface  layer  will  be  modeled  by  solving 
the  corresponding  Navier-Stokes  equations  for  steady- 
state  linear  flow  explicitly.  The  subsequent  analysis  leads  to 
an  explicit  expression  for  the  drilling  velocity.  The  drilling 
depth  is  then  calculated  by  numerical  integration  of  the 
drilling  velocity  and  compared  to  published  experimental 
data.  Finally,  the  relationship  between  drilling  depth  and 
recoil  momentum  is  discussed. 


II.  EVAPORATION 

The  rapid  surface  evaporation  due  to  absorption  of 
intense  laser  radiation  has  been  modeled  by  various 
authors.*""  It  occurs  whenever  the  absorbed  energy  den¬ 
sity  is  much  greater  than  the  amount  necessary  to  evapo¬ 
rate  the  target  material. 

The  description  of  Knight®  shall  be  used  in  this  work, 
since  it  considers  evaporation  into  an  outside  atmosphere 
rather  than  into  vacuum.  The  inclusion  of  the  outside  at¬ 
mosphere  leads  to  an  enhancement  of  the  forces  induced  by 
the  evaporated  material  onto  the  target  surface,  as  one 
might  intuitively  guess.  Knight’s  description  assumes  a 
quasiequilibrium  evaporation  from  a  plane  liquid-vapor  in¬ 
terface,  the  evajwration  front,  into  air  at  ambient  condi¬ 
tions.  The  model  considers  a  Knudsen  layer  extending 
some  molecular  mean  free  paths  from  the  liquid  into  the 
vapor  region."  It  assumes  a  discontinuous  change  of  tem¬ 
perature  r,  pressure  p,  and  density  p  across  the  evapora¬ 
tion  front.  The  flow  of  the  expanding  hot  vapor,  which 
drives  a  shock  wave  into  the  surrounding  air,  determines 
the  strength  of  these  discontinuities.  They  are  given  in  Ref 
8  as 
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where  the  subscript  v  denotes  the  variables  describing  the 
vapor  adjacent  to  the  Knudsen  layer  and  s  the  saturated 
water  vapor.  y=  Cp/c„  stands  for  the  ratio  of  heat  capaci¬ 
ties.  denotes  the  flow  Mach  number  of  the  vapor  leav¬ 
ing  the  Knudsen  layer.  Sonic  flow  applies  for  the  case  of 
evaporation  into  vacuum.^  For  evaporation  into  an  outside 
atmosphere,  the  flow  is  sonic  only  if  the  incident  laser 
fluxes  are  high  enough.  Then  equals  1,  while  for  sub¬ 
sonic  flow  it  follows  from  Ref.  8  as 


Pa 


6lll 


(4) 


Here,  pg,  and  denote  specific  beat  ratio,  ambient 
pressure,  and  sound  speed  in  air;  y^=1.4  and 
=  where  /?„=  ^  is  the  universal  gas 

constant,  and  is  the  molecular  weight;  =  29 
g/gmol  for  air.  Further,  the  sound  speed  in  water  vapor, 
«»=  is  calculated  using  y=  j  and  18 

g/gmol. 

Equations  (l)-(4)  contain  eight  unknowns;  pressure 
p,  temperature  T,  and  density  p,  for  saturated  and  expand¬ 
ing  vapor,  plus  the  Mach  flow  numbers  m  and  Mg.  With 
the  aid  of  the  equation  of  state,  p—pRT,  for  saturated  and 
expanding  vapor,  these  variables  can  all  be  expressed  as  a 
function  of  the  saturated  vapor  pressure  p,  and  tempera¬ 
ture  Tr  These  two  variables  are  related  by  a  Clausius- 
Clapeyron  relationship  as 

(5) 


A  and  T.  were  determined  by  fitting  published  vapor  pres¬ 
sure  data  of  water  with  Eq.  (S)  (Ref.  13);  using  the  values 
A  =  3.0705 X 10^  bars  and  r.  =  4688.2  K,  Eq.  (5)  fits  the 
data  in  the  range  firom  120  *C  to  the  critical  temperature 
Tg  =  374  *C,  with  an  accuracy  better  than  2%. 

Energy  balance  at  the  evaporation  front  provides  the 
missing  equation  to  determine  the  thermodynamic  state  of 
the  saturated  vapor.  To  formulate  it  I  introduce  the  vari¬ 
able  the  pait  of  the  laser  irradiance  that  is 

actually  used  to  evaporate  target  material.  The  relationship 
between  7*,^  and  will  be  discussed  below.  The  location 
of  the  liquid-vapor  interface  is  denoted  by  D(t).  The  ve¬ 
locity  with  which  it  moves  into  the  target  as  a  consequence 
of  evaporative  mass  removal  will  be  referred  to  as  the  ther¬ 
mal  drilling  velocity  Dg  Using  this  notation,  the  energy 
balance  for  a  small  volume  fixed  at  the  evaporation  front, 
at  z  =  D(t),  reads 


+  Pii>tfil=PgM^g^Hg  +  ^  , 

where  the  specific  enthalpy  H  =€  +  p/p,  e  being  the  inter¬ 
nal  energy,  has  been  introduced.  Note  that  energy  is  added 
to  the  small  volume  not  only  by  the  laser  beam,  but  also  by 
the  transport  of  the  internal  energy  of  the  liquid  into  it, 
represented  by  the  second  term  on  the  left  side  of  Eq.  (6). 
Using  mass  conservation  pfi,=  pMfiv  across  the  evapora¬ 
tion  front,  the  energy  balance  is  rewritten  as 

Iey*p=Pfifllu<  (2) 

where  the  notation  has  been  further  simplified  by  introduc¬ 
tion  of  7f/„(7’„)=  (7f„  —  Hi)  +  yA7j/?u,r,/2.  In  this  sum 
the  enthalpy  diflerence  //„  —  Hi  is  much  larger  than  the 
kinetic-energy  term.  For  an  initial  temperature  of  20  °C, 
one  calculates,  using  published  data,*^ 
Hg  -  77,  =  2.630  ±  0.088  kJ/g  and  =  2.73  ±0.1kJ/g  in 
the  temperature  range  from  90  to  340  °C. 

To  connect  Eq.  (7)  with  the  saturated  vapor  pressure 
Pp  I  write  the  equation  of  state  of  the  saturated  vapor  in  the 
form 

A = ZiT  P»  ^  ylr^wT,  J  ^  y/rJ^u,Tg.  (8) 

yMg  pg  Mg 

Using  mass  balance  across  the  evaporation  front  to¬ 
gether  with  the  energy  balance  equation  (7),  this  gives 

(9) 

The  force  per  area  po  that  acts  on  the  surface  of  the 
liquid  is  obtained  from  momentum  conservation  across  the 
evaporation  front.  It  is  equal  to  the  momentum  flux  of 
vapor  from  the  Knudsen  layer,  p,(pg/pt)  (1-1-  yMl).  Us¬ 
ing  the  equation  of  state  of  the  vapor  and  Eq.  (9)  gives 

(10) 

Note  that  here  the  momentum  of  the  liquid  relative  to 
the  moving  evaporation  front,  pb^„  has  been  neglected.  The 
force  acting  on  the  surface  thus  becomes  a  product  of  the 
evaporated  mass  per  unit  time,  the  sound  speed  in  satu¬ 
rated  vapor,  and  the  expression  ylTg/Tg{M\ 
-K  \/y)/Mg  that  is  dominated  by  its  dependence  on  the 
Mach  number  Mg.  For  water  at  an  initial  temperature  of 
20 ’C  and  a  1-atm  back  pressure,  sonic  flow  applies  for 
surface  temperatures  7',>223*C.  Numerically,  one  then 
finds  V77r,(A7j  -|-  l/y)/M,.=  1  578. 

For  high  enough  irradiances,  7ev»p>10^  W/cm^,  the 
pressure  Po  is  proportional  to  It  is  therefore  useful  to 
define  the  proportionality  factor  as 

J‘rv^^evap')~Po^^cvap'  (11) 

It  is  plotted  in  Fig.  1(a)  as  a  function  of  /evap»  P®ri  of 
the  laser  irradiance  that  is  used  for  evaporation.  The  vari¬ 
ation  of  is  mainly  caused  by  the  changing  Mach  number 
Mg  and  is  pronounced  only  at  irradiances  smaller  than  lo’ 
W/cm^. 
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(a)  Irradiance,  I„,p  [MW/cm^] 


FIG.  1.  Descriptioii  of  pure  evaporation  as  a  function  of  the  irradiance 
spent  on  evaporation,  /^.  The  target  is  assumed  to  be  water  at  norma! 
conditions.  The  onset  of  sonic  flow  at  =  0.44  MW/cm^  is  indicated 
by  dashed  lines,  (a)  Ratio  of  pressure  acting  on  the  target  surface  and 
irradiance  spent  on  evaporation,  /„=  vs  The  drop  of  the 

curve  at  low  irradiances  is  mainly  due  to  the  changing  Mach  number,  (b) 
Saturated  vapor  temperature  T,  vs  T,  varies  very  roughly  logarith¬ 
mically  with  irradiance.  Sonic  flow  applies  for  r,>223  'C 

Because  of  the  exponential  dependence  in  Eq.  (S),  the 
temperature  T,  depends  approximately  logarithmically  on 
the  irradiance  This  is  plotted  in  Fig.  1(b)  for  tem¬ 
peratures  below  the  critical  temperature  of  water, 
r,  =  374*C 

Hi.  EJECTION 

The  preceding  analysis  of  evaporation  assumed  a  con¬ 
stant  irradiance  beam  profile.  In  what  follows,  a  TEMoq 
laser  beam  is  assumed.  To  describe  this  I  apply  the  results 
obtained  above  locally.  This  is  a  reasonable  approximation 
for  describing  the  drilling  velocity,  because  the  Knudsen 
layer  is  very  thin  compared  to  the' focal  spot  size,  and  for 
the  case  of  sonic  flow,  which  occurs  at  high  enough  inci¬ 
dent  irradiances,  the  determination  of  the  surface  condi¬ 
tions  is  independent  of  the  exterior  vapor  flow. 

As  the  pressure  po  increases  monotonically  with  the 
irradiance  spatial  variations  of  the  Gaussian  beam 


Laser 

Im 

i  I  I  I  I  I  M 


vapor 


FIG.  2.  Illustration  of  the  drilling  mechanism.  Between  vapor  and  solid  a 
liquid  layer  of  thickness  6  forms.  The  induced  material  flow  enhanres 
drilling.  At  the  center  of  the  crater,  radial  flow  is  a  good  approximatioD. 


intensity  generate  radial  pressure  gradients  acting  on  the 
irradiated  surface.  If  the  shear  forces  generated  within  that 
surface  layer  are  stronger  than  the  shear  strength  of  the 
material  in  that  layer,  it  will  be  displaced  radially  (cf.  Fig. 
2).  This  can  be  expected  if  vaporization  of  the  target  ma¬ 
terial  is  preceded  by  liquefaction,  as  is  assumed  here.  The 
resulting  liquid  flow  is  determined  by  the  evaporation-in¬ 
duced  radial  pressure  gradient  as  well  as  by  the  thickness 
and  viscosity  of  the  liquid  layer.  At  the  bottom  of  the 
crater,  the  material  is  displaced  predominantly  in  radial 
direction.  The  crater  walls  deflect  the  flow  into  axial  direc¬ 
tion.  This  eventually  leads  to  the  observed  expulsion  of 
liquid  from  the  drilled  holes.  The  generation  of  a  liquid 
layer  thus  enhances — and  complicates — the  ablation  pro¬ 
cess  considerably. 

A.  The  thickness  of  the  liquid  layer 

The  thickness  of  the  liquid  layer  that  separates  solid 
from  gas  at  the  evaporation  front  depends  on  the  depth  of 
energy  deposition.  For  ablation  of  gelatin  by  infrared  radi¬ 
ation,  this  is  determined  by  the  optical  penetration  depth 
rather  than  by  heat  conduction.  Indeed,  the  typical  dis¬ 
tance  be  over  which  conduction  acts  is  yfiftc, 
K  =  1.44X  10  “  ^  cmVs  being  the  heat  diffusivity  for  water. 
The  time  during  which  conduction  acts  equals  roughly 
the  time  it  takes  to  ablate  a  layer  of  thickness  bg,  i.e.,  the 
thickness  bg  divided  by  the  drilling  velocity  D.  Thus,  if  the 
layer  thickness  during  ablation  was  due  to  heat  conduction 
alone,  6c  =  k/D,  less  than  about  0. 1  ^m  for  drilling  veloc¬ 
ities  larger  than  1  m/s.  Compared  with  the  minimal  optical 
penetration  depth  of  water  in  the  IR,  d  =  0,75  ^m,’*'**  this 
is  indeed  negligible.  Note  that  for  metals  this  is  not  so;  for 
liquid  aluminum  and  a  drilling  velocity  of  1  m/s,  one  cal¬ 
culates  K/b  =  39  nm,  which  is  much  more  than  the  optical 
penetration  depth. 
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Temperature,  Tj 


FIG.  3.  Illustration  of  the  temperature  distribution  Tj  caused  by  irradia¬ 
tion  of  the  surface  of  the  target  at  ^  =  0  from  above  [cf.  Eq.  (13)].  Close 
to  the  surface  T,  is  claniped  at  T,  the  saturated  vapor  temperature. 
Deeper  inside  the  target,  it  drops  off  exponentially  with  the  distance 
from  the  surface. 


The  extension  of  the  liquid  layer  can  be  calculated 
approximately  by  applying  energy  balance  to  an  infinites¬ 
imal  volume  element  at  a  fixed  distance  ^  =  z  —  D(t)  away 
from  the  evaporation  front.  In  case  of  pure  evaporation,  no 
radial  material  flow  exists.  Then,  the  energy  balance  near 
the  solid-liquid  boundary  reads 


where  p/=  1  g/cm^  and  c,  =  4.18  J/(g  K)  are  the  density 
and  specific  heat  of  the  liquid,  the  values  for  water  being 
assumed. 

Assuming  the  radiation  to  decay  exponentially 
with  a  penetration  depth  d,  the  temperature  of  the  liquid, 
Tf,  is  calculated  for  an  initial  temperature  7,  =  lO’C  as 

T^^)~Ti=T,-T,  for 


^enp 

t: — ^ 

Dp/Ci 


for 


(13) 


The  resulting  temperature  distribution  is  plotted  in  Fig.  3. 
For  temperatures  below  the  saturated  vapor  temperature 
Tp  T I  decreases  exponentially  with  the  distance  from  the 
surface,  f.  Ti  is  equal  to  the  melting  temperature  T„  at  a 
distance  away  from  the  surface  and  to  T„  at  another 
distance  In  this  description  the  liquid  temperature  Tj 
cannot  rise  higher  than  at  distances  closer  to  the  sur¬ 
face,  The  excess  absorbed  energy  in  that  region  is 
assumed  to  be  spent  entirely  on  the  liquid-vapor  phase 
transition. 

With  the  aid  of  Eq.  (7),  I  obtain  from  the  above  Eq. 
(13): 


U=d\n 


H,, 

c,(T„-T,)- 


(14) 


This  expression  for  the  liquid-layer  thickness  was  derived 
assuming  purely  evaporative  drilling.  In  the  case  to  be  con¬ 
sidered,  where  ejection  of  liquid  increases  the  ablation 
speed,  the  actual  temperature  decay  will  be  steeper  and  the 
resulting  thickness  of  the  liquid  layer  somewhat  smaller 
than  the  prediction  of  Eq.  (14).  However,  in  view  of 
the  logarithmic  dependence,  the  associated  changes  are 
small. 

The  liquid  flow  also  depends  strongly  on  the  kinematic 
viscosity  of  the  molten  target  material,  which  for  gelatin  is 
difficult  to  access  experimentally  above  100  "C,  as  at  these 
temperatures  it  decomposes  within  a  few  minutes.  Around 
100  °C  I  have  measured  the  viscosity  and  found  that  it 
varies  about  exponentially  with  temperature,  as 
v(7)=voexp[  —  v/(7  —  ■&)].  I  shall  apply  this  approx¬ 
imation  to  extrapolate  the  viscosity  of  gelatin  to  higher 
temperatures.  To  calculate  the  drilling  depth  in  gelatin 
containing  0.84  g/cm^  water  numerically,  I  shall  use  below 
Vo  =  0.11  cmVs,  V,  =  0.0513 ‘C“',  and  0  =  97.5‘'C  (cf. 
Sec.  IV). 

To  simplify  the  subsequent  analysis  of  the  flow,  I  shall 
assume  that  the  temperature  inside  the  liquid  layer  is  equal 
to  some  average  value  7,v=  7,  -f  a,v(7,  —  7,).  The  pa¬ 
rameter  compares  the  average  rise  of  temperature  of 
the  displaced  liquid  to  the  one  at  the  surface  of  the  target. 
It  depends  on  the  axial  profiles  of  temperature  and  radial 
flow  velocity  within  the  liquid.  Here,  it  will  be  assumed 
that  a.v  is  constant. 

Correspondingly,  the  thickness  8  of  the  liquid  layer  is 
approximated  as  (cf.  Fig.  3) 


6=dln 


Hi, 

cPtATs-  Ti)  ■ 


(15) 


For  the  experimental  situation  in  mind,  the  Uquid  layer  8 
has  a  thickness  of  about  1.5  optical  penetration  depths  or 
1.1  fim.  Note  that  in  this  approximation  the  viscosity 
within  the  liquid  layer  becomes  equal  to  an  average  value, 
independent  of  f . 


B.  The  ejection  mechanism 

The  laser-induced  pressure  gradients  at  the  evapora¬ 
tion  front  push  material  radially  outward,  away  from  the 
axis  of  the  driUing  beam.  The  displaced  material  can  then 
either  flow  along  the  walls  and  out  of  the  crater  or,  as 
observed  with  targets  that  deform  during  ablation,  remain 
trapped  within.^  Either  way,  as  long  as  the  displaced  ma¬ 
terial  does  not  obstruct  the  central  part  of  the  beam,  the 
drilling  velocity  at  the  center  of  the  drilled  hole  increases. 
The  total  velocity  of  the  evaporation  front  into  the  target, 
D,  thus  becomes  the  sum  of  two  contributions:  a  thermal 
drilling  velocity  caused  by  evaporation.  Dp  and  a  mechan¬ 
ical  drilling  component  caused  by  material  displacement, 
£)„.  Its  dependence  on  the  irradiance  of  the  drilling  beam 
will  now  be  discussed. 

To  describe  the  liquid  flow,  I  assume  that  the  liquid 
layer  is  planar  and  its  thickness  small  compared  to  the 
beam  waist  (Oq  of  the  drilling  beam.  The  Navier-Stokes 
equations  then  take  the  form  of  boundary-layer  equations; 
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dVr  dv^  dVr  1  dp 
dp 

0  =  ^.  (17) 

where  v,  and  v^  denote  the  radial  and  axial  flow  velocities. 
The  continuity  equation  for  cylindrical  flow  is 

dv.  1  dv. 

Dimensional  analysis  shows  that  the  two  nonlinear  terms 
in  the  Navier-Stokes  equation  (16)  have  the  same  order  of 
magnitude.  They  can  be  neglected  if  v^v/dr-<  vd^v/d^, 
which  is  correct  for  slow  enough  flow,  i.e.,  u^<<v/6. 

Steady-state  flow  determines  the  resulting  hole  depth, 
provided  that  it  develops  within  a  time  i,  that  is  short 
compared  to  the  laser-pulse  duration.  For  pulse  durations 
longer  than  100  /xs,  this  is  the  case,  as  can  be  seen  from 
dimensional  considerations.  Indeed,  if  one  assumes  that  all 
the  contributions  in  Eq.  (16)  are  the  same  order  of  mag¬ 
nitude,  comparing  the  acceleration  term  with  the  viscosity 
term  gives  y,/tj=  or  t,a-6Vv.  Introducing  the  ap¬ 

propriate  numbers  for  gelatin  at  100  *C,  t,«0.1  fis. 

I  therefore  proceed  to  solve  the  linearized  steady-state 
Navier-Stokes  equations 

1  dp  d^Vr 


As  these  equations  are  linear,  one  can  easily  solve  for  the 
flow  velocities  in  a  moving  coordinate  system,  setring 
^  =  z  —  D{t)  (compare  Fig.  2).  At  the  liquid-solid  inter¬ 
face  f  =  6,  no  radial  flow  is  allowed,  implying 
u,(^  =  6)  =  0.  At  f  =  0,  the  free  surface  of  the  liquid,  the 
correct  boundary  condition  is 

P=Po,  (21) 

dv, 

v^=0,  (22) 

where  pg  is  the  evaporation-induced  pressure,  given  by  Eq. 
(10).  The  axial  flow  velocities  at  the  phase  interfaces  do 
not  vanish.  In  the  moving  frame  of  reference,  they  are 
i;,(i:  =  0)=:  -  i),  and  i;,(^  =  6)  =  -  (i>,  +  £>„)■ 

The  solution  to  this  flow  problem  is 

I  dpo  .  , 

<23) 

I  can  now  give  an  expression  for  b„,  the  mechanical 
contribution  to  drilling.  Integrating  the  equation  of  conti¬ 
nuity  (18)  gives 


lim 

Va-  ._o 


1  d \dpQ 
r  ^  dr)  dr 


Dimensional  analysis  of  Eq.  (25)  shows  that  the  mechan¬ 
ical  drilling  velocity  increases  with  growing  layer  thickness 
S,  decreasing  viscosity  v,  or  increasing  radial  surface  pres¬ 
sure  gradient  roughly  as  b„c^p^^/(p(V(o\).  This  result 
differs  from  Eq.  (15)  in  Ref.  6  by  the  dimensionless  factor 
8/ciQ.  The  difference  arises  because  in  Eq.  (12)  of  that 
paper  the  expression  used  for  the  frictional  force  is  errone¬ 
ous.'^  Nevertheless,  the  general  conclusions  drawn  there 
remain  valid. 

Equation  (25)  remains  accurate  as  long  as  the  linear 
approximation  holds,  i.e.,  v^Kv/S.  As  the  axial  flow  veloc¬ 
ity  is  of  the  same  order  of  magnitude  as  the  drilling 
velocity,  the  linear  approximation  is  accurate  for  small 
drilling  velocities  only.  For  Ji)>v/6,  Eq.  (25)  is  but  an 
upper  limit  to  the  mechanical  drilling  velocity.  The  non¬ 
linear  inertia  terms  in  the  Navier-Stokes  equations  (16) 
and  (17)  lead  to  a  smaller  steady-state  mechanical  drilling 
velocity  b„  as  they  add  resistance  to  the  flow.  It  can  be 
shown  that  for  highly  nonlinear  flow  the  steady-state  drill 
ing  velocity  is  given  by  b„  i=:4{8/ci)q)  yjpo/pt-  This  descrip¬ 
tion  of  the  drilling  mechanism  then  becomes  equivalent  to 
the  one  given  by  von  Allmen*  that  was  based  on  thermo¬ 
dynamic  arguments  alone.  It  is  interesting  to  note  that,  as 
Pq  is  about  proportional  to  the  mechanical  drilling 
velocity  b„  scales  differently  wi»h  /evap.  the  irradiance  used 
for  evaporation,  in  the  two  liquid-flow  regimes. 

Independent  of  the  details  of  the  flow,  the  mechanical 
contribution  to  drilling  increases  as  the  spot  size  is  re¬ 
duced.  Note,  however,  that  this  increase  is  limited  by  en¬ 
ergy  conservation,  because  as  the  liquid  flow  carries  more 
heat  out  of  the  irradiated  region,  less  energy  will  be  avail¬ 
able  for  the  liquid-vapor  phase  transition.  As  a  conse¬ 
quence  evaporation,  the  driving  mechanism  of  ejection, 
lessens  and  the  surface  temperature  T,  drops. 

The  drilling  velocity  b  can  be  calculated  with  the  aid 
of  Eqs.  (7),  (11),  and  (25)  once  the  dependence  of  the 
irradiance  consumed  by  evaporation,  /ev»p.  on  the  incident 
irradiance  has  been  established.  To  this  end  I  consider 
the  energy  balance  for  a  small  volume  centered  on  axis: 

+  P/^'/«»v(2's  -  Ti)b„.  (26) 

Because  b„  depends  on  derivatives  of  the  pressure  pg,  Eq. 
(26)  is  actually  a  differential  equation.  To  obtain  an  ap¬ 
proximate  solution,  I  assume  /ev,p(r)  to  be  proportional  to 
Iin('’).  the  incident  irradiance.  Tliis  assumption  is,  obvi¬ 
ously,  the  more  accurate  the  smaller  b„.  Now  for  a 
Gaussian  beam,  the  mechanical  drilling  velocity  is  calcu¬ 
lated  from  Eqs.  (11)  and  ( 25 )  as 


8  6  -' 


■^cvap 


/cv  +  /, 


'"•p  dl„ 


—  f"  +  r  -  (24)  and /ev».-('' =  0)  follows  numerically  from  Eq.  (26).  Note 


Using  the  above  boundary  conditions  and  Eq.  (23),  this 
integration  is  easy.  It  follows  that 


that  for  lcnf>0. 1  MW/cm^  the  proportionality  factor  /„ 
is  about  constant  and,  therefore,  the  mechanical  drilling 
velocity  proportional  to  pg  [cf.  Fig.  1(a)]. 
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It  is  useful,  at  this  point,  to  simplify  the  notation  by 
introduction  of  the  drilling  velocity  ratio  q)  and  the  energy 
consumption  parameter  ip  as 

<P=bjb,  (28) 


(29) 


The  drilling  velocity  ratio  q>  describes  the  hydrodynamic 
aspect  of  ablation.  For  linear  flow  it  follows  from  Eqs,  (7) 
and  (25).  Its  variation  with  is  mainly  caused  by  the 
strong  temperature  dependence  of  the  viscosity  and,  for 
subsonic  gas  flow,  by  the  influence  of  the  Mach  number  on 
the  proportionality  factor  The  consumption  parameter 
Ip  compares  the  cost  of  the  two  drilling  processes  in  terms 
of  energy.  It  depends  but  weakly  on  the  irradiation  condi¬ 
tions. 

Using  this  notation,  I  solve  Eqs.  (7)  and  (26)  for  b,  to 
obtain  the  total  drilling  velocity 


1  I. 

D=- 


1  -I-  <p 


p,H,„  \-\~(pip' 


(30) 


As  ^<1,  the  drilling  velocity  grows  monotonically  with 
the  drilling  velocity  ratic  q>.  If  ejection  becomes  dominant, 
^>1,  the  drilling  velocity  b  approaches  the  value  /;„/ 
[p</a,v(7’,  —  r,)].  This  still  depends  implicitly  on  tp,  be¬ 
cause  the  amount  of  energy  spent  on  the  liquid-vapor  phase 
transition  determines  the  temperature  T,  even  if  that  en¬ 
ergy  is  negligible  when  compared  to  the  heat  carried  away 
by  the  flowing  liquid. 


IV.  DISCUSSION 


In  this  section  the  model  will  be  used  to  discuss  drilling 
depth  and  recoil  momentum,  two  physical  parameters  that 
i  yield  important  information  on  the  drilling  mechanism. 

They  both  depend  on  the  initial  position  Dq  of  the  target 
i  surface  with  respect  to  the  beam  waist  in  a  characteristic 

[way  that  is  qualitatively  very  similar  for  metals  and  bio¬ 
logical  targets.''^’* 

Published  measurements,  obtained  from  ablation  of 
gelatin  containing  0.84  g/cm^  water  with  single 
;YAG  laser  pulses  at  a  wavelength  of  A  =  2.94  ^m, 
will  now  be  discussed  in  detail.  In  these  experiments  a 
TEMoo  beam  with  a  beam  waist  of  <wo  =  72  (im  and  a  peak 
irradiance  of  2.45  MW/cm^  was  used.  The  laser  pulses  had 
a  duration  of  t  =  200  /xs  and  a  pulse  energy  of  £„  =  40  mJ. 
These  parameter  values  will  be  assumed  in  the  following. 
Note  that  the  peak  fluences  used  in  these  experiments, 
~490  J/cm^,  are  much  larger  than  the  threshold  value  for 
,  ablation,  which  for  gelatin  is  similar  to  guinea  pig  skin,  ~  1 
J/cm^” 


\  A.  brilling  depth 


If  the  experimentally  determined  drilling  depth  is 
small  compared  to  the  confocal  length  Zo=  irox^/X  of  the 
focused  laser  beam,  the  irradiation  conditions  at  the  evap¬ 
oration  front  do  not  change  during  ablation.  The  resulting 
drilling  depth  Z),  then  just  equals  the  drilling  velocity  b 
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FIG.  4.  Influence  of  the  distance  between  the  surface  of  the  target  and  the 
oeara  waist  oi.  the  drilling  mechanism.  Positive  (negative)  distance  values 
correspond  to  a  bea.ti  waisl  position  located  outside  (inside)tbe  target. 
Experimental  data  from  drilling  gelatin  targets  containing  0.84  g/cm^ 
water  with  single  laser  pulses  is  denoted  by  solid  circles  (•)  (from  Ref. 
1 ).  The  corresponding  laser  parameters  were  A  =  2.94  /rm,  t  =  200  /is, 
<U(,  =  72  /im,  and  =  40  mJ.  (a)  Drilling  depth  vs  waist-surface  dis¬ 
tance.  Drilling  depth  predicted  by  the  model  using  <z„  =  0.75  (solid) .  For 
comparison  the  calculated  drilling  depth  in  case  of  pure  evaporation  is 
also  shown  (dashed),  (b)  Recoil  momentum  vs  waist-surface  distance. 
The  line  is  to  guide  the  eye.  The  relative  minimum  of  the  recoil  is  ob¬ 
served  at  the  same  initial  waist-surface  distance  as  the  maximum  of  the 
drilling  depth.  This  is  a  consequence  of  minimizing  the  amount  of  evap¬ 
orated  mass. 


multiplied  by  the  pulse  duration  t.  Under  tight  focusing 
conditions,  however,  it  is  quite  common  to  drill  to  depths 
comparable  with  or  larger  than  the  confocal  length.  In  that 
case  the  drilling  velocity  changes  with  the  irradiance  at  the 
moving  evaporation  front,  and  the  hole  depth  D^.  must  be 
calculated  by  numerical  integration.  Equivalently,  one  can 
solve  numerically  the  integral  equation 


rAj-t  A 
Jbo 


(31) 


for  Dr  In  Fig.  4(a)  the  dependence  of  the  total  drilling 
depth  on  the  initial  position  Dg  is  shown.  The  solid  line 
is  the  result  of  a  calculation  that  uses  Eq.  (30)  to  express 
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the  drilling  velocity  and  =  0.75.  This  value  for  the 
only  remaining  free  parameter  in  the  calculation,  was  ob¬ 
tained  by  fitting  the  maximal  calculated  depth  to  the  cor¬ 
responding  experimental  value.  It  follows  that  the  average 
temperature  increase  of  the  displaced  liquid  at  the  center  of 
the  crater  is  about  75%  of  the  surface  temperature  rise 
T  —T 

Also  shown  in  Fig.  4(a)  is  the  calculated  drilling  depth 
versus  target  position  curve  for  pure  evaporation  (^  =  0), 
given  by  the  dashed  line.  One  notes  that  the  maximum 
drilling  depth  is  more  than  3.5  times  larger  than  that  which 
would  result  from  pure  evaporation.  As  the  target  is  shifted 
far  away  from  the  maximum  drilling  position,  ejection 
ceases  to  be  important.  Only  then  can  the  experimental 
drilling  data  be  explained  by  purely  evaporative  drilling. 

The  model  overestimates  the  drilling  depth  in  the 
wings  of  the  curve.  The  parameter  a.v  affects  mainly  the 
maximum  drilling  depth;  changing  it  does  not  improve  the 
fit.  Indeed,  the  width  of  the  calculated  curve  in  essence 
depends  on  the  product  8^/{p(vo?)  [cf.  Eq.  (25)].  There¬ 
fore,  the  deviation  of  the  calculated  depths  from  the  exper¬ 
imental  results  seems  to  be  caused  by  the  large  uncertainty 
in  the  viscosity  of  molten  gelatin  at  temperatures  substan¬ 
tially  above  100  ”C  and  the  omission  of  the  depth  depen¬ 
dence  of  the  viscosity  within  the  liquid  layer.  The  reason 
why  this  inaccuracy  does  not  affect  the  calculation  more 
seriously  is  that  the  drilling  velocity  is  controlled  to  a  large 
extent  by  energy  conservation  at  the  evaporation  front. 
Moreover,  the  omission  of  the  nonlinear  terms  in  the  de¬ 
scription  of  the  liquid  flow  as  well  as  the  one-dimensional 
treatment  of  the  vapor  flow  should  be  kept  in  mind. 

The  dependence  of  drilling  velocity  on  position  z  is 
symmetric  with  respect  to  the  beam  waist  position,  where 
it  attains  a  maximal  value.  On  the  other  hand,  the  drilling 
depth  versus  initial  target  position  curve  cannot  be  sym¬ 
metric  because  the  ablation  velocity  rises  as  one  drills  to¬ 
ward  the  waist  while  it  lessens,  once  that  position  is  passed. 
Note  also  that  the  model  predicts  quite  high  temperatures 
in  the  center  of  the  crater  for  the  experimental  conditions 
used,  7’,=:200‘’C. 

B.  Recoil  momentum 

The  recoil  momentum  transferred  to  the  target  during 
ablation  is  often  used  to  characterize  laser-material  inter¬ 
actions,  as  it  is  easily  accessible  to  measurement,  e.g.,  in  a 
ballistic  pendulum  arrangement.  It  is  usually  described  by 
the  mechanical  coupling  coefficient  C,  which  is  defined  as 
the  ratio  of  the  transferred  momentum  to  the  delivered 
pulse  energy  and  has  units  of  dyn  s/J.  For  ablation  of  alu¬ 
minum  with  laser  pulses  of  a  duration  of  1  /xs,  the  mechan¬ 
ical  coupling  has  been  shown  to  result  from  bulk  target 
vaporization.^® 

In  Fig.  4(b)  the  experimentally  determined  variation 
of  the  coupling  coefficient  C  with  the  distance  between  the 
target  surface  and  the  beam  waist  is  shown.  It  is  very  in¬ 
teresting  to  relate  the  experimental  recoil  and  drilling- 
depth  data  to  each  other.  A  coarse  estimate  of  the  recoil  is 
obtained  by  assuming  that  all  the  ablated  material  leaves 
the  target  perpendicular  to  its  surface,  the  evaporated  ma¬ 


terial  at  a  speed  equal  to  and  the  liquid  ejecta  at  an 
average  “ejection  velocity”  Ui,q,  much  less  than  the  speed  of 
sound.  The  time  derivadve  of  the  coupling  coefficient  can 
thus  be  estimated  by  the  sum  of  the  respective  products  of 
speed  and  mass  removed  per  time  as 

dC  Pprcol  . 

(32) 

Using  Eqs.  (28)  and  (30),  this  is  rewritten  as 

For  large  waist-surface  distances,  the  irradiated  spot  is 
big  and  ejection  negligible,  ^  =  0.  Using  Eq.  (33)  and  as¬ 
suming  sonic  flow,  the  calculated  coupling  coefficient  then 
becomes  C  ~  aJHi^  At  a  vapor  temperature  T„oi  150*C, 
this  is  much  higher,  C~  19  dyn  s/J,  than  what  is  actually 
observed.^' 

The  discrepancy  arises  because  the  vapor  flow  was 
modeled  one  dimensionally,  implying  a  flow  perpendicular 
to  the  target  surface.  However,  molecules  evaporated  from 
different  regions  of  the  irradiated  spot  can  interact  with 
each  other  as  soon  as  they  have  traveled  a  distance  com¬ 
parable  to  the  spot  diameter.  As  they  leave  the  target  at 
about  the  speed  of  sound,  this  happens  after  a  very  short 
fraction  of  the  pulse  duration  used,  t  =  200  /xs.  Therefore, 
to  accurately  describe  the  recoil  caused  by  the  vapor,  its 
flow  must  be  modeled  two  dimensionally.  Furthermore,  it 
has  been  shown  that  the  presence  of  the  surrounding  air  as 
well  as  the  finite  surface  area  of  the  irradiated  target  influ¬ 
ence  the  mechanical  coupling  coefficient.^^’^^ 

Still,  for  large-aspect-ratio  holes,  the  gas  flow  is  quite 
directional  independent  of  pulse  duration  and  the  model 
here  presented  helps  in  qualitatively  understanding  the  ex¬ 
perimental  data.  Inspection  of  Eq.  (33)  shows  that  in  that 
case  the  time  derivative  of  the  coupling  coefficient  de¬ 
creases  as  the  velocity  ratio  grows,  only  if 

<\f>,  being  in  the  range  from  0.1  to  0.3  for  tem¬ 
peratures  T,  between  100  and  300  "C.  This  is  a  reasonable 
assumption,  because  the  liquid  slows  down  as  it  flows  along 
the  crater  wall.  Also,  some  of  the  displaced  liquid  remains 
in  the  crater,  reducing  the  average  ejection  velocity  Uijq 
further. 

For  aluminum  the  loss  of  mass  of  the  targets  has  been 
determined  simultaneously  with  the  coupling  coefficient  as 
a  function  of  the  distance  between  the  target  surface  and 
the  beam  waist.  In  these  experiments  the  loss  of  mass  was 
maximal,  like  the  drilling  depth,  at  the  waist-target  dis¬ 
tance  that  corresponded  to  the  minimal  coupling  coeffi¬ 
cient.^  When  drilling  in  gelatin,  the  loss  of  mass  and  the 
drilling  depth  do  not  correlate  as  simply;  under  certain 
irradiation  conditions  a  large  part  of  the  displaced  liquid 
remains  trapped  inside  the  crater  because  of  elastic  defor¬ 
mations  that  occur  during  ablation.’ 

Assuming  U]jq/(Af,xi„)  <^,  Eq.  (33)  thus  implies  that 
as  ejection  is  enhanced,  the  coupling  coefficient  should  de¬ 
crease  because  of  the  reduction  of  the  amount  of  evapo¬ 
rated  mass.  For  the  long  pulse  durations  in  mind,  the  in¬ 
fluence  of  aspect  ratio  and  total  mass  evaporated  on  the 
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target  recoil  can  therefore  be  interpreted  as  follows:  For 
shallow  holes  the  vapor  flows  mainly  parallel  to  the  surface 
of  the  target,  with  but  a  small  perpendicular  velocity  com¬ 
ponent.  As  the  radial  momentum  components  cancel  each 
other,  the  coupling  coeflScient  C  becomes  considerably 
smaller  than  that  which  would  be  predicted  by  Eq.  (33). 
Increasing  the  aspect  ratio  raises  the  average  velocity  com¬ 
ponent  of  the  vapor  flow  perpendicular  to  the  surface  of  the 
target.  Consequently,  the  coupling  coefiicient  increases,  as 
observed  in  Fig.  4(b)  for  target  positions  corresponding  to 
drilling  depths  less  than  about  1  mm.  Drilling  deeper  holes 
leads  to  smaller  coupling  coefficients.  This  is  because  the 
further  increase  of  the  hole  depth  reduces  the  total  mass 
evaporated,  but  does  not  enhance  the  directionality  of  the 
vapor  flow  substantially.  The  coupling  coefficient  thus  de¬ 
creases. 

The  minimum  in  the  curve  shown  in  Fig.  4(b)  is  thus 
a  consequence  of  minimizing  the  total  amount  of  mass 
evaporated  during  the  ablation  process.  This  is  equivalent, 
however,  to  maximizing  the  mechanical  drilling.  As  dis¬ 
placing  material  is  more  efficient  with  respect  to  drilling 
depth  than  evaporating  it,  it  becomes  clear,  therefore  why 
the  maximal  drilling  depth  and  the  relative  minimum  of 
the  coupling  coefficient  occur  at  the  same  target  position. 

V.  CONCLUSION 

Liquefaction  of  surface  material  prior  to  evaporation 
changes  the  laser-target  interaction  qualitatively.  It  en¬ 
hances  drilling  efficiency  and  leads  to  a  reduction  of  the 
evaporated  mass  and,  consequently,  of  the  surface  temper¬ 
ature.  These  effects  are  important  for  ablation  using  pulses 
of  a  duration  of  100  /rs  or  longer  under  tight-focusing  con¬ 
ditions.  They  become  insignificant  if  the  liquid  cannot 
transverse  the  irradiated  spot  during  the  laser  pulse,  as  is 
typically  the  case  for  ^-switched  pulses. 

A  steady-state  model  is  developed  that  describes  abla¬ 
tion  as  a  combination  of  evaporation  and  ejection.  Using 
energy  balance  and  linearized  Navier-Stokes  equations  for 
a  boundary  layer,  an  upper  limiting  value  for  the  drilling 
velocity  is  obtained.  The  drilling  velocity  is  found,  as  a 
rule,  to  increase  with  increasing  depth  of  energy  deposi¬ 
tion,  decreasing  liquid  viscosity,  and  smaller  spot  size. 

The  recoil  momentum  induced  on  the  target  depends 
predominantly  on  the  vapor-flow  pattern.  This  becomes 


more  and  more  directional  as  the  aspect  ratio  of  the  drilled 
holes  increases.  In  that  case  the  coupling  coefficient  is  dom¬ 
inated  by  the  total  amount  of  evaporated  mass. 
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Abstract 

Ablation  of  biological  targets  by  highly-absorbed  in¬ 
frared  laser  pulses  (duraticm  200  fts)  is  expired 
using  a  model  that  takes  into  account  tissue  Uque- 
faction.  The  generation  of  a  liquid  layer  in  tissue 
leads  to  an  enhancement  of  the  drilling  efficiency  and 
to  a  reduction  of  the  evaporation  rate.  Experimen¬ 
tal  drilling-depth  and  recoil-momentum  data  are  ex¬ 
plained  by  a  combination  of  evaporation  and  ejection. 


1  Introduction 

In  modem  surgery  COj  lasers  are  being  used  rou¬ 
tinely  to  perform  high  predsion  cutting  of  tissue.  To¬ 
day  the  Er:  YAG  laser,  which  emits  at  the  infrared  ab¬ 
sorption  peak  of  water,  allows,  in  principle,  an  even 
higher  cutting  precision.  However,  with  the  appli¬ 
cation  of  pulsed  lasers  in  surgery  there  has  been  a 
growing  awareness  that  tissue  ablation  is  an  "explo- 
ave  event”  rather  than  a  straightforward  evaporative 
process  [1]. 

There  has  been  no  direct  experimental  evidence 
for  the  existence  of  a  liquid  zone  during  ablation. 
This  is  primarily  because  it  is  thin,  of  the  order  of 
the  optical  penetration  depth  of  the  laser  radiation, 
and  hidden  inside  the  drilled  crater.  However,  in¬ 
direct  confirmation  of  the  liquid  layer  is  provided 
by  some  experimental  findings  obtained  using  laser 
pulses  longer  than  100  fis  under  tight-focusing  condi¬ 
tions,  i.e.  spot  diameters  of  less  than  half  a  mm.  Of 
particular  practical  importance  are  the  extension  and 
repetition-rate  dependence  of  lateral  thermal  damage 
along  COj  laser  incisions  as  well  as  the  axial  trans¬ 
port  of  biologically  active  material  into  laser  craters 
12]  [3]  [4]. 

Here  1  intend  to  further  solidify  the  concept  of  “li¬ 
quefied  tissue”  by  connecting  recent  drilling-depth 
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and  recoil-momentum  measurements  on  tissue  and 
gelatin  targets.  After  a  description  of  the  generation 
of  the  liquid  layer  1  proceed  to  explain  how  ablation 
can  be  understood  as  a  combination  of  evaporation 
and  hquid  ejection.  This  concept  is  then  applied  to 
drilling  depth  and  recoil  momentum. 


2  Liquid  Layer 

As  a  target  absorbs  intense  radiation,  a  surface  layer 
is  heated  and  begins  to  evaporate.  Even  before  eva¬ 
poration  that  layer  undergoes  structural  changes.  Its 
thickness  depends  on  the  depth  of  energy  deposi¬ 
tion.  For  infrared-laser  tissue  ablation  the  thickness  is 
usually  determined  by  absorption  of  radiation  rather 
than  heat  conduction. 

Suppose  that  the  thickness  of  the  liquid  layer  was 
determined  predominatly  by  heat  conduction.  The 
distance,  £«,  over  which  conduction  acts  is  y/i^, 
where  k  «  1.22  •  10~®  cm^/s  is  the  thermal  ffiffu- 
sivity  for  excised  animal  dermis  [5].  After  the  on¬ 
set  of  drilling,  the  time,  tc,  during  which  conduc¬ 
tion  acts  is  given  by  the  time  it  takes  to  ablate  a 
layer  of  thickness,  6^,  rather  than  simply  the  laser 
pulse  duration.  In  that  case  tc  becomes  equal  to  the 
layer  thickness,  6c,  divided  by  the  drilling  velocity,  D. 
The  layer  thickness  caused  by  heat  conduction  alone 
thus  amounts  to  =  k/D,  which  is  already  equal 
to  the  minimal  optical  penetration  depth  in  the  in¬ 
frared,  d  =  0.75  fiin  at  2.94  /rm  wavelength  [6]  [7], 
at  a  drilling  velocity  of  0.16  m/s.  For  larger  drUling 
velocities,  therefore,  the  extension  of  the  liquid  layer 
is  dominated  by  the  radiation  absorption. 

The  extension  of  the  liquid  layer  can  be  calculated 
by  applying  energy  balance  to  an  infinitesimal  volume 
element  at  a  fixed  distance  away  from  the  surface  of 
the  liquid  layer.  Assuming  a  temperature  of  the  liquid 
of  150  ®C  one  finds  that  the  thickness  of  the  liquid 
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Figure  1:  Photomicrograph  showing  the  ultrastructure  of  collagen  fibrils  in  the  coriura  of  (ICR)  mice  in  cross 
(C)  and  longitudinal  (L)  section.  Freshly  excised  skin  samples  were  thermally  altered  by  iiiunersion  in  a  heat 
bath  of  temperature  7*  for  160  seconds,  (a):  T*  =  53  *C:  The  fibrils  are  normal,  (b):  Tj  =  65  ®C:  The  fibrils 
are  thickened  to  about  twice  the  normal  size,  (c)  Tj  =  80  ®C:  The  collagen  fibres  are  completely  dissolved. 
X  31. 000  (from  [3]) 


layer  generated  during  ablation  is  about  1.5  times  the 
optical  penetration  depth. 

Water  and  gelatin  melt  at  well  defined  tempera¬ 
tures  which  can  easily  be  measured.  It  is  less  obvi¬ 
ous  how  a  liquid  layer  is  generated  in  tissue,  which 
does  not  undergo  an  equally  clear-cut  phase  transi¬ 
tion.  Tissue  can  be  regarded  as  “liquefied”,  however, 
as  soon  as  its  shear  strength  is  overcome  by  shear 
forces.  Then  adjacent  layers  of  tissue  move  relative 
to  each  other,  much  as  in  any  viscous  flow.  For  in¬ 
tact  skin  this  occurs  at  pressures  comparable  to  its 
tensile  strength,  100  bar  [8].  Such  high  ablation 
pressures  are  not  generated  at  irradiances  of  less  than 
~  1  MW/cm^. 

Thermally  damaged  tissue  is  expected  to  flow  at 
much  lower  pressures,  as  the  associated  structural 
changes  weaken  the  tensile  strength  considerably.  In 
skin,  for  example,  where  collagen  fibers  are  respon¬ 
sible  for  the  high  tensile  strength,  this  assumption 
is  supported  by  an  observation  of  their  thermally- 
induced  decay.  In  Fig.  1  the  ultrastructure  of  colla¬ 
gen  fibrils  harvested  from  thermally  altered  mice  skin 
(corium)  is  shown.  The  skin  samples  were  thermally 
damaged  by  immersion  in  a  heat  bath  of  a  constant 
temperature  for  160  seconds.  Fig.  1  shows  how  the 
collagen  fibrils  dissolve,  and  therefore  weaken,  with 
increasing  temperature,  Tj,  of  the  heat  bath.  Long 
(160  seconds)  exposure  at  80  *C  destroys  the  collagen 
structure  completely.  Similar  changes  are  expected 
for  short  thermal  exposures,  provided  that  the  tem¬ 


peratures  are  high  enough. 

In  summary,  during  ablation  there  is  a  surface  layer 
that  has  a  reduced  tensile  strength  and  flows  more 
easily  under  the  action  of  shear  forces  than  the  bulk  of 
the  target.  Its  thickness  is  of  the  order  of  the  optical 
penetration  depth. 


3  Ablation  Model 

Theoretical  models  of  ablation  tend  to  become  elabo¬ 
rate  when  pushed  to  increasing  rigorousness.  I  will 
describe  the  ablation  process  qualitatively,  stressing 
the  principles  involved  rather  than  their  mathemati¬ 
cal  description.  A  more  detailed  presentation  of  simi¬ 
lar  ideas  can  be  found  in  the  literature  [9]. 

Vaporization  of  materials  by  intense  laser  radiation 
has  been  modelled  by  many.  Early  works  treat  va¬ 
porization  into  vacuum  [10]  [11],  while  a  more  recent 
one  also  considers  the  influence  of  a  non-vanishing 
back  pressure  [12].  One  main  conclusion  of  these 
works  is  that  the  surface  temperature  increases  mono- 
tonically  but  weakly  with  the  incident  irradiance. 
This  b  found,  in  principle,  by  combining  a  Clausius- 
Qapciron  equation  for  saturated  vapor  temperature 
and  pressure  with  the  fact  that  the  vapor  pressure  de¬ 
pends  linearly  on  the  amount  of  evaporated  material 
(cf.  also  [13]). 

The  pressure  that  pushes  on  the  surface  from  which 
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Figme  2:  Schematic  <liagiam  illufitrating  the  model. 
At  the  bottom  of  the  crater,  which  has  a  depth,  D(t), 
a  liquid  layer  of  thickness,  6,  forms  between  vapor 
and  solid.  Material  is  pushed  outwards,  enhancing 
drilling.  At  the  center  of  the  crater  the  flow  is  to  a 
good  approximation  radial. 


evaporation  proceeds  is  determined  by  momentum 
conservation.  It  increases  with  irradiance  and  is 
about  equal  to  55%  of  the  vapor  pressure  for  sonic 
flow  [14].  The  induced  ablation  pressures  will  there¬ 
fore  change  spatially  with  the  variation  of  irradiance 
on  the  target.  This  induces  shear  forces  inside  the  li¬ 
quid  surface  layer.  Thus  for  a  Gaussian  beam,  where 
the  intensity  drops  from  the  center  to  the  edge,  li¬ 
quid  flows  radially  towards  the  periphery  (cf.  Fig.  2). 
Other  beam  profiles  may  cause  a  more  complicated 
flow  structure  within  the  liquid  layer.  However,  some 
Uquid  will  always  be  pushed  towards  the  crater  wall. 
At  the  wall  the  liquid  changes  direction  and  flows 
axially  out  of  the  hole.  How  much  of  this  displaced 
material  actually  leaves  the  crater  depends  on  its  flow 
properties  as  well  as  on  the  shape  of  the  crater  wall 
during  ablation.  Such  questions  are  important  with 
respect  to  the  thermal  damage  remaining  after  abla¬ 
tion,  but  will  not  be  addressed  here. 

To  quantify  the  enhancement  of  drilling  by  this 
mechanical  ejection  mechanism,  one  needs  to  study 
the  hydrodynamics  of  the  induced  flow.  Here  I  shall 
take,  for  simplicity,  a  rather  phenomenological  point 
of  view  and  assume  that  ejection  is  present,  and  that 
its  effect  is  determined  from  experiment.  The  drilling 
velocity,  D,  thus  becomes  a  sum  of  two  terms,  the 
thermal  drilling  velocity  caused  by  evaporation,  £>«, 
and  the  mechanical  contribution.  Dm-  The  strength 
of  the  mechanical  drilling  component  is  conveniently 


measured  in  terms  of  the  corresponding  evaporative 
contribution.  Therefore  I  define  the  ratio  of  the  two 
drilling  velocities  (at  the  center  of  the  crater) 

<p=bmlDt  (1) 

<p  depends  on  the  irradiation  parameters  as  well  as  on 
the  hydrodynamic  properties  of  the  target  material. 
Note  th*U,  <p  does  not  directly  give  the  amplification  of 
the  drilling  velocity  due  to  the  presence  of  the  liquid 
layer,  but  rather  the  relative  strength  of  the  two  con¬ 
tributions,  under  the  assumption  that  they  are  both 
present.  <p  increases  as  it  becomes  easier  to  move 
liquid  in  the  radial  direction.  Therefore,  for  given 
irradiance,  ip  grows  with  increasing  thickness  and  de¬ 
creasing  viscosity  of  the  liquid  layer,  as  well  as  with 
decreasing  spot  size  of  the  drilling  beam.  This  is  only 
true  for  sufficiently  long  pulses.  Otherwise  the  liquid 
moves  less  than  a  spot  radius,  wq,  during  the  pulse 
duration,  and  the  nnechanical  contribution.  Dm,  be¬ 
comes  negligible.  This  will  be  the  case  for  Q-switched 
pulses  with  a  duration  of  r  10  ns.  Taking  the  sound 
speed  in  water,  a  =  1.5  km/s,  as  an  upper  limit  for 
the  flow  velocity,  the  product,  ra  «  15  pm,  is  smaller 
than  most  spot  sizes  realized  in  ablation  experiments. 

A  good  qualitative  understanding  of  the  ablation 
process  can  be  gained  from  an  energy  balance  at 
the  surface  of  the  irradiated  target.  I  assume  steady 
state  drilling  has  been  established.  The  incident  ir¬ 
radiation,  I  in,  is  then  split  between  the  two  energy¬ 
consuming  processes,  evaporation  and  ejection. 

lin  =  piHivbt  -1-  piCiiJ'i  —  Ti)bm  (2) 

Here  pt,  Cf  and  Ti  denote  the  density,  the  specific 
heat  and  the  temperature  of  the  liquid.  I  shall  use 
thermal  properties  of  water,  i.e.  pj  =  1  g/cm*,  cj  = 
4.2  J/(g  ■  K)  and  an  initial  temperature  of  7)  =  20  "C. 
Hi„  is  the  enthalpy  difference  between  the  target  ma¬ 
terial  at  the  initial  temperature,  7),  and  the  evapo¬ 
rated  target  material.  Hi„  depends  only  weakly  on 
the  temperature  of  the  vapor;  =  2.73  ±  0.1  kJ/g 
in  the  temperature  range  from  100  **0  to  the  critical 
temperature,  T*  =  374  ®C. 

Using  Eqs.  (1-2)  the  total  drilling  velocity  D  = 
Dt  -f  Dm  can  be  written  as 

b=  ^-ti¬ 
pi  ffiv  +  •pci{Ti  —  Ti) 


According  to  equation  (3)  the  drilling  velocity  in¬ 
creases  steadily  with  the  contribution  of  ejection. 

=  0  characterizes  the  purely  evaporative  process, 
while  =  oo  gives  the  upper  limit  on  the  drilling 
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Velocity  Ratio,  cp 


Figure  3:  Calculated  liquid  temperature,  7),  ver¬ 
sus  (p  —  Dm/Dt,  the  ratio  between  the  mechanical 
drilling  velocity,  £>„,  and  the  evaporative  contribu¬ 
tion,  Dt,  for  an  incident  irradiance  of  2.5  MW/cm^. 
The  temperature,  Tj,  is  determined  predominantly 
by  the  evaporation  rate.  Therefore  it  drops  with 
increasing  <p,  as  more  and  more  energy  is  used  to 
heat  up  liquid  displaced  during  the  ablation.  At 
<p  ~  2.5  (dashed  line),  the  two  drilling  processes  con¬ 
sume  equal  amounts  of  energy. 

velocity,  when  the  energy  spent  on  evaporation  can 
be  neglected.  If  y?,  the  ratio  between  mechanical 
and  thermal  drilling  velocity,  increases,  less  energy 
is  available  for  evaporation.  The  evaporation  rate 
is  thus  reduced;  this  lowers  the  liquid  temperature, 
T|.  This  is  illustrated  in  Fig.  3.  There  the  liquid 
temperature,  T|,  is  plotted  versus  fp.  The  calculation 
assumed  a  target  irradiance  of  2.5  MW/cm^.  For  a 
small  mechanical  attribution  the  energy  is  mainly 
used  to  evaporate  material,  while  othervrise  a  signifi¬ 
cant  amount  of  energy  is  spent  to  heat  up  the  liquid 
displaced  during  ablation.  The  breakeven  point,  de¬ 
termined  from  Fig.  3  and  eq.  (2),  is  found  at  y?  =  2.5. 

4  Discussion 

In  this  section  the  model  will  be  used  to  explain 
drilling  depth  and  recoil  measurements.  Both  are 
quite  easily  accessible  to  erq)erimental  investigation 
and  yield  important  information  on  the  drilling  mech¬ 
anism.  Published  measurements,  obtained  from  abla¬ 
tion  of  tissue  and  gelatin  targets  with  single  Er ;  YAG 
laser  pulses  will  be  discussed  [15].  In  these  exper¬ 
iments  a  TEMoo-beam  with  a  beam  wrust  of  u^o  = 
72  ftm  and  a  peak  irradiance  of  2.5  MW/cm*  was 


used.  The  laser  pulses  had  a  duration  of  r  =  200  fis 
and  a  pulse  energy  of  £■,„  =  40  mJ.  These  parameter 
values  will  be  assumed  in  the  following. 

4.1  Drilling  depth 

Under  tight-focusing  conditions  it  is  quite  common  to 
drill  to  depths  comparable  to  or  larger  than  the  confo- 
cal  distance,  zo  =  xwq/A.  Then  the  irradiance  at  the 
phase  boundary,  and  therefore  also  the  drilling  velo¬ 
city,  change  significantly  during  drilling.  As  a  conse¬ 
quence  the  total  drilling  depth  is  not  just  a  product 
of  the  pulse  duration  and  the  initial  surface  recession 
speed,  but  has  to  be  calculated,  instead,  by  numerical 
integration  using  the  drilling  velocity. 

The  drilling  depth  becomes  maximal  if  the  beam 
waist  is  located  exactly  halfway  down  the  drilled 
crater.  A  little  thought  shows  that  this  holds  be¬ 
cause  of  the  monotonic  increase  of  the  drilling  velo¬ 
city  with  irradiance.  In  Fig.  4  the  calculated  maximal 
drilling  depth  is  plotted  versus  the  ratio,  ip,  of  the 
two  material  removing  processes  assuming  the  laser 
parameters  given  above.  The  dashed  lines  indicate 
the  maximal  depths  in  skin  and  gelatin,  determined 
from  experimental  data  plotted  in  Fig.  5a,  as  well  as 
the  corresponding  average  values  of  <p.  This  calcula¬ 
tion  is  approximate,  as  it  assumes  <p  is  constant,  while 
in  reality  ip  depends  on  irradiance  as  well  as  on  the 
details  of  the  flow  and  changes  during  the  ablation. 
Still  the  calculation  illustrates  the  consequence  of  the 
ejection  mechanism  on  the  maximal  drilling  depth. 

In  Fig.  5a  the  drilling  depth  measured  in  tissue  and 
in  gelatin  containing  80%  water  is  displayed  as  a  func¬ 
tion  of  the  distance  between  the  beam  waist  and  the 
target  surface  at  the  beginning  of  the  laser  pulse.  The 
sho.rt-dashed  line  indicates  the  (calculated)  depth  for 
purely  evaporative  drilling  into  a  water  target.  As 
is  made  clear  by  the  flgure,  the  experimental  values 
are  too  large  to  be  caused  by  evaporation  only.  The 
enhanced  drilling  depth  originates  from  the  ejection 
mechanism. 

4.2  Recoil  momentum 

The  recoil  momentum,  caused  by  the  ejected  mate¬ 
rial  during  an  ablation  process  on  a  target,  provides 
another  simple  means  of  characterizing  the  interac¬ 
tion,  as  it  can  be  measured  quite  easily  in  a  ballistic 
pendulum  arrangement.  The  total  recoil  momentum 
divided  by  the  pulse  energy  is  called  the  mechani¬ 
cal  coupling  coefficient,  C,  and  has  units  of  dyn  •  s/J. 
In  Fig.  5b  its  dependence  on  the  beam  waist  to  tar- 
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Figure  4:  Calculated  maximal  drilling  depth  ver¬ 
sus  ip  =  DmiDt,  the  ratio  between  the  mechanical 
drilling  velocity,  and  the  evaporative  contribu¬ 
tion,  Ai  for  an  incident  irradiance  of  2.5  MW/cm^. 
Depths  were  calculated  by  numerical  integration  of 
drilling  velocity  assuming  the  beam  waist  to  be 
halfway  down  the  laser  crater.  The  dashed  lines  indi¬ 
cate  the  experimentally  determined  maximal  depths 
in  skin  and  gelatin  containing  0.84  g/cm^  water  and 
the  corresponding  average  <p  values  (cf.  Fig.  5a). 

get  distance  is  shown.  The  similar  appearance  of  the 
curves  for  gelatin  and  skin  indicates  that  the  laser- 
target  interaction  is  qualitatively  identical  for  the  two 
media. 

To  roughly  estimate  the  total  recoil  momentum 
caused  by  the  ablation  process,  I  assume  that  ejection 
contributes  little.  Indeed  this  seems  to  be  a  reason¬ 
able  assumption,  as  the  liquid  is  slowed  down  during 
its  flow  along  the  crater  walls.  It  will  therefore  leave 
the  target  at  a  much  lower  speed  than  the  vapor.  Fur¬ 
thermore,  not  all  of  the  displaced  liquid  is  actually 
ejected  during  the  ablation. 

Therefore,  neglecting  the  contribution  of  liquid 
ejecta,  the  recoil  momentum  is  estimated  by  assuming 
that  all  the  evaporated  material  leaves  the  target  at 
the  speed  of  sound,  a,.  Correspondingly  the  coupling 
coefficient,  C,  becomes 

C  «  (p/  rbt  )a,  lEin  (4) 


Figure  5:  Influence  of  the  initial  position  of  the  target 
surface  on  the  drilling  mechanisni.  Positive  (negative) 
positions  correspond  to  a  beam  waist  located  outside 
(inside)  the  target.  Elxporimental  data  is  indicated  by 
open  circles  (o)  for  skin  targets  and  by  filled  circles  (•) 
for  gelatin  containing  0.84  g/cm*  water.  All  the  holes 
were  drilled  with  single  laser  pulses  (A  =  2.94  pm, 
T  =  200  ps.  Wo  =  72  pm,  Ein  =  40  mJ).  Data 
from  [15].  The  hues  are  intended  to  guide  the  eye. 
(a):  Drilhng  depth  versus  waist-target  distance.  For 
comparison  the  drilling  depth  in  case  of  pure  evapo¬ 
ration,  calculated  by  numerical  integration,  is  shown 
as  a  short-dashed  line,  (b):  Coupling  coefficient,  C, 
versus  waist-target  distance.  C  is  equal  to  the  total 
recoil  momentum  normalized  by  the  laser  pulse  en¬ 
ergy.  Note  that  C  has  a  relative  minimum  when  the 
drilhng  depth  has  a  maximum. 


Using  Eqs.  (1,3)  this  is  rewritten  as 

H,„+<pc,{Ti-Ti) 

For  large  spot  sizes,  where  there  is  no  ejection,  e.g. 
(p  =  0,  the  calculated  coupling  coefficient  becomes 


C  r;  OvlHiv  At  a  vapor  temperature  of  Tp  =  150  ®C 
this  is  much  higher,  C  «  19  dyn  •  s/3,  than  what  is 
actually  observed,  <  3  dyn  s/J  *. 

’  Taking  into  account  the  back  pressure  of  the  outside  at- 
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The  discrepco'-y  arises  because  the  vapor  flow  was 
modelled  one-  dimensionally,  implying  a  flow  perpen¬ 
dicular  to  the  target  surface.  However,  molecules 
evaporated  from  different  regions  of  the  irradiated 
spot  cam  interact  with  each  other  as  soon  as  they  have 
travelled  a  distance  comparable  to  the  spot  diameter. 
As  they  leave  the  target  at  about  the  speed  of  sound, 
this  happens  already  after  a  very  short  fraction  of  the 
pulse  durations  used,  r  =  200  fis.  Therefore,  to  ac¬ 
curately  describe  the  recoil  caused  by  the  vapor,  its 
flow  must  be  modelled  two-dimensionally  [16].  Nev¬ 
ertheless,  for  large  aspect-ratio  holes  the  gas  flow  is 
quite  directional  independent  of  pulse  duration.  Then 
the  model  here  presented  helps  in  qualitatively  un¬ 
derstanding  the  experimental  data.  It  predicts  that 
the  total  amount  of  mass  evaporated  determines  the 
recoil  momentum. 

The  recoil  curves  can  be  interpreted  by  the  in¬ 
fluence  of  aspect  ratio  and  total  mass  evaporated: 
For  shallow  holes  the  vapor  flows  mainly  parallel  to 
the  tauget  surface,  with  only  a  small  velocity  com¬ 
ponent  perpendicular  to  the  surface  of  the  target. 
As  the  radial  momentum  components  cancel  each 
other  because  of  symmetry,  the  coupling  coefficient, 
C,  thus  becomes  considerably  smaller  than  the  pre¬ 
diction  based  on  the  eq.  (5).  Increasing  the  aspect 
ratio  raises  the  perpendicular  velocity  component  of 
the  vapor  flow,  and  therefore  also  the  coupling  coeffi¬ 
cient.  This  rise  is  seen  in  Fig.  5b  for  target  positions 
corresponding  to  drilling  depths  smaller  than  about 
a  mm.  For  larger  drilling  depths  the  coupling  coeffi¬ 
cient  is  seen  to  decrease  with  increasing  depths.  This 
is  because  the  further  increase  of  the  hole  depth  does 
not  enhance  the  directionality  of  the  vapor  flow  sub¬ 
stantially,  but  it  reduces  the  total  mass  evaporated. 
The  coupling  coefficient  thus  decreases,  as  predicted 
by  eq.  (5)  due  to  the  larger  average  tp  value  associ¬ 
ated  with  the  enhanced  mechanical  contribution  to 
drilling. 

The  minima  in  the  curves  shown  in  Fig.  5b  are  thus 
a  consequence  of  minimizing  the  total  amount  of  mass 
evaporated  during  the  ablation  process.  This  is  equi¬ 
valent,  however,  to  maximizing  the  ejection  mecha¬ 
nism.  As  ejection  is  more  efficient  than  evaporation 
with  respect  to  drilling  depth,  it  becomes  clear,  there¬ 
fore,  why  the  maximal  drilling  depth  and  the  relative 
minimum  of  the  coupling  coefficient  in  Fig.  5  occur 
at  the  same  target  position. 


■noephere  does  not  help  enough,  as  the  associated  reduction  of 
the  coupling  coefficient,  — Peii//  «  —1  dyn  •  s/J,  is  too  small. 


5  Conclusion 

Liquefaction  of  the  target  surface  prior  to  evapo¬ 
ration  changes  the  laser-target  interaction  qualita¬ 
tively.  The  associated  effects  are  important  for  ab¬ 
lation  of  tissue  by  pulses  of  a  duration  of  100  fts  or 
longer  under  tight-focusing  conditions.  They  become 
insignificant  if  the  liquid  cannot  transverse  the  ra¬ 
dius  of  the  irradiated  spot  during  the  laser  pulse,  as 
is  typically  the  case  for  Q-switched  pulses. 

Drilling  depth  and  recoil  momentum  caused  by 
pulsed  laser  ablation  change  with  distance  of  the  tar¬ 
get  surface  from  the  beam  waist  position  in  a  cha¬ 
racteristic  way,  qualitatively  identical  for  tissue  and 
gelatin  targets.  The  experimental  findings  are  ex¬ 
plained  by  a  model  that  describes  ablation  as  a  com¬ 
bination  of  evaporation  and  ejection. 

The  generation  of  a  liquid  layer  at  the  phase  boun¬ 
dary  and  its  radial  displacement  during  ablation  lead 
to  an  increase  of  the  drilling  efficiency.  The  associated 
reduction  of  the  evaporation  rate  causes  a  reduction 
of  the  temperature  at  the  target  surface  during  abla¬ 
tion. 

The  recoil  momentum  induced  on  the  target  de¬ 
pends  predominantly  on  the  vapor  flow  pattern.  This 
becomes  more  and  more  directional  as  the  aspect  ra¬ 
tio  of  the  drilled  holes  increases.  In  that  case  the 
coupling  coefficient  is  dominated  by  the  total  amount 
of  evaporated  mass. 
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sure  wave)  while  the  lower  trace  displays 
the  Ho:YAG  pulse  waveform.  The  firet 
pressure  pulse  corresponds  to  the  jjeak  of 
the  laser  pulse  and  is  believed  to  result 
fiom  the  mechanism  of  opttical  sound  gen¬ 
eration.'  The  second  and  higher-magni¬ 
tude  pressure  pulse  occurs  at  the  end  of 
the  Ho;YAG  ktser  pulse,  coinciding  with 
the  timing  for  bubble  collapse  (Rg.  1) 
which  is  believed  to  be  the  source  for  this 
pressure  wave.  The  subsequent  pressure 
signals  may  be  caused  by  the  ringing  of 
the  transducer  or  by  echoes,  and  needs 
further  clarification. 

Generation  of  cavitation  bubbles  and 
the  associated  pressure  transients  may  be 
important  mechanisms  for  the  tissue  tear¬ 
ing  aitd  dissection  observed  from  previ¬ 
ously  reported  mid-infrared  User  soft  tis¬ 
sue  ablation  studies.'  In  addition,  the 
strong  interaction  between  mid-infrared 
laser  pulses  with  fluid  and  the  induced 
pressure  wave  should  be  carefully  moni¬ 
tored  in  ’ophthalnuc  use  of  these  lasers. 
Understandme  and  cf>nrml  nt  -res- 
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CThL4  Shock  waves  generated  by  XeQ 
excimer  laser  abUtion  of  polyimide  in 
air  and  water 

A.  D.  Zweig,  T.  F.  Deutsch,  Wellman 
Leboratories  of  Photomedicine,  Massachusetts 
General  Hospital,  Boston,  MA  02U4 

Pressure  transients  of  several  hundred 
bars  with  rise  times  of  a  few  nanoseconds 
can  be  created  by  UV-User  abUtion  of 
organic  polymers  as  well  as  of  tissue*'' 
We  demonstrate,  we  believe  for  the  first 
time,  that  ablation  of  polyimide  by 
30B-nm  XeO  excimer  laser  radiation  gen- 
oates  shock  waves  at  fluenoes  of  several 
times  the  threshold  of  abUtion.  This  is  a 
significant  fiitding,  because  a  shock  wave, 
characterized  by  a  picosecond  rise  time 
and  a  supersoitic  propagation  speed,  dif¬ 
fers  qualitatively  from  an  acotistic  wave. 
Our  measurements  were  intended  to  sim- 
uUte  the  abUtion  of  cornea  at  193  nm,  but 
are  also  relevant  to  material  processing  of 
lasers. 

We  investigated  the  abUtion  of  sheet 
poiyimide  fkMting  on  water  with  a  free 
sutUkx  in  air,  as  well  as  of  polyimide 
inunersed  in  distilled  water.  The  stress 
waves  that  are  Uunched  in  the  polyintide 
during  ^^tion  propagate  through  the 
25iim  thick  sam|te  aitd,  after  a  r^ection 
loss,  into  dte  water  at  the  rear  side  of  the 
target.  There  we  determirte  their  propa¬ 
gation  speed  from  a  measurement  of  the 
time  delay  between  the  deflection  of  two 
parallel,  focused,  helium-neon  probe 
beams  that  ate  separated  by  —40  pm. 
The  vdodty  measurement  suffices  to 
diaracterize  tlw  shock,  as  the  leUtionship 
between  shock  velocitv  and  pressure  is 


At  fluences  of  only  a  few  times  the 
threshold  of  abUtion  F„  the  ablabon-in- 
duced  pressure  transients  propagate  at 
supersonic  velocities  in  the  water,  indica¬ 
ting  the  generation  of  shock  waves.  The 
peak  pressures  of  the  shocks  propagating 
in  the  water  are  measured  near  the  rear 
surface  of  the  target.  Pressure  P. 
increases  with  fluence  F  as  (F— PJ’,  where 
n  =  2/3  if  the  abUtion  is  in  air  and  n  = 
1/2  if  the  abUtion  is  confined  by  water. 
Figures  1  and  2  show  the  experimental 
results. 

Also  shown  in  these  figures  are  the  the¬ 
oretical  predictions  of  an  ideal  gas  model 
that  was  used  earlier  by  other  workers  to 
explain  shock  generation  by  confined 
abUtion  of  metals.*'*  According  to  that 
model  the  pressure  depends  only  on  irra- 
dUnce  and  acoustic  impedances  of  the 
materials  involved,  provided  that  the 
optical  absorption  coefficient  is  high 
enough.  It  is  seen  that  the  mcxlel 
describes  the  pressures  generated  by  the 
confined  abUtion  quite  successfully. 

For  abUtion  in  air  (Hg.  1)  the  pressures 
predicted  by  the  model  are  lower  than  the 
experiment^  results.  This  may  be  due  to 
subsurface  heating;  if  subsurface  layers 
are  heated  sufficiently  rapidly,  their 
expansion  can  be  confined  by  the  overly¬ 
ing  material.  This  implies  that  even  the 
abUtion  of  a  free  surface  may  actually  be 
confined  significantly  by  the  target  mate¬ 
rial  itself.  The  confinement  effect  would 
increase  the  abUtion  pressures.  If  the 
energy  were  deposited  deeper  in  the 
polyimide  sample,  possibly  by  using 
longer  wavelen^  radUtion,  the  abUtion 
pressures  could  in  principle  rise  even 
more.  The  corresponding  model  predic¬ 
tion  is  indicated  in  Hg.  1  as  a  dash^  line. 

Although  these  measurements  were 
designed  to  measure  the  pressures  accom¬ 
panying  pulsed  excimer  abUtion,  it  is 
worth  noting  that  the  same  abUtion  sys¬ 
tem  may  now  be  used  as  a  source  of  char¬ 
acterized  shock  waves.  Such  a  source 
may  have  use  in  sonochemistry,  as  well  as 
in  the  study  of  shexk  effects  on  cells  and 
macromole^es. 

1.  P.  E.  Dyer  and  R.  Srinivasan.  Appl 
Phys.  Lett  48,  445  (1986). 
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Dyer.  J.  Appl,  Phys.  64,  2194  (1988) 
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CThLS  Time  resolved  transmission  of 
collagen-based  films  daring  ArF  excimer 
laser  ablation 

M.  N.  Ediger,  G.  H.  Pettit,  R.  P.  Weiblingor, 
FDA  Center  for  Devices  &  Radiological 
Health,  1901  Chapman  Ave.,  Rockville,  MD 
20857 

Time  resolved  reflectivity  studies  of  ArF 
(193-tun)  laser  pulses  on  bovine  cornea 
have  shown  substantial  temr>oral  trunca 


results  have  also  been  obtained  using  cii 
Ugen  gel  substrates.'  These  investig,- 
tions  have  indicated  that  this  pulse"'  ..j, 
tion  is  primarily  due  to  dynamic  n. 
ity  changes  in  the  siuface,  and  is  nc 
predominantly  an  effect  of  scattering  b 
abUtion  products.  Since  the  underlyin; 
physical  principles  of  pulsed,  ultraviolc 
User  tissue  ablation  are  not  yet  wel 
understood,  dynamic  changes  in  the  opti 
cal  properties  of  the  cornea  are  of  interes 
in  the  study  of  the  energy  deposition  am 
tirartsport.  In  this  work,  the  time  resolve^ 
transmission  of  oolUgen  filnts  and  of  thii 
sections  of  bovine  cornea  have  been  mea 
sured  to  provide  mformation  about  th( 
light/tissue  interaction  at  the  surface  o 
collagen-based  substrates. 

Thin  filrtxs  (<10  jun)  of  denatured  colU 
gen  (geUtin)  were  irradUted  with  the  spa 
tially  filtered  outpmt  of  an  ArF  excimei 
User.  A  typrical  transmitted  wavefom 
occuring  at  a  User  fluence  of  110  m]  •  cm~ 
is  comprared  with  the  incident  laser  pmlsc 
in  Rg.  1.  In  this  figure,  both  waveform; 
have  been  normalized  to  equal  p>eak  mag 
nitude.  Clearly,  the  transmission  of  the 
sample  increases  over  the  course  of  tht 
irradiation  pjulse.  This  is  in  distinci 
oppxrsition  to  subablative  irradUtior 
where  the  trarrsmitted  waveform  exhibiti 
exact  correspxrndence  to  the  pulse  shapx 
of  the  incident  User  waveform.  Artalysu 
is  presented  that  indicates  that  the  tempx> 
ral  distortion  at  abUtive  fluences  is  no' 
simply  due  to  a  reduction  in  the  thicknes; 
of  the  sample  by  abUtion.  The  a<  aa 
tion  of  the  trailing  edge  of  the  tiarv  ..tec 
pulse  correUtes  with  reflectivity  measure 
ments  in  colUgen  gels.*  In  that  work,  th( 
reflecti'vity  was  observed  to  fall  precipi 
tously  in  the  tail  of  the  pulse,  suggestinp 
energy  deposition  and/or  transmission  o 
the  sample  were  radically  increasing  a 
the  same  irtsfant. 

The  results  of  transmission  expreriment 
prerformed  with  thin  sections  (—10  jim)  o 
bovine  cornea  are  shown  in  Hg.  2 
Again,  the  mddent  User  waveform  i; 
shown  for  comprarison,  and  the  pnilse 
have  been  nonnalized  to  ecpual  magni 
tilde.  The  User  fluence  for  corneal  filn 
irradUtion  was  170  mj  •  cm''.  The  bans 
nutted  waveform  for  the  comeal  film  dif 
fers  considerably  in  shapre  from  tha 
obtained  for  geUtin  films.  Instead,  th< 
pulse  transnvitted  by  comeal  films  mor 
resembles  the  truncated  reflected  pniU 
observed  from  whole,  unsectioited  cor 
neas.  These  results  suggest  that  th' 
attenuation  of  the  User  by  the  intact  bf 
sue  irKteases  over  the  course  of  the  abla 
tion  pmlse. 

The  dissirttilarity  in  the  transmitted  pml 
ses  of  these  two  oolUgen-based  substrate 
may  be  due  to  the  confirmational  diffei 
ences  in  the  absorbing  protein.  TSes 
expjerimental  results  are  prresented  an 
contrasted  with  those  obtairred  usin 
films  prroduced  with  bovine  coUager  m  ii 
native  (nondenatiired)  confirmatioi 

1.  G  H.  Pettit,  M.N.  Ediger,  and  R. 
Weiblinger,  “Dynamic  opjtical  propx' 
ties  of  colUgen-based  tissue  durin 
ArF  excimer  l^r  abUtion",  submittc 
to  Appl.  Opt.,  Nov.  1991 . 

2.  M.  N.  Ediger,  G.  H.  Pettit,  and  R. 
Weiblinger,  “"Bovine  comeal  phon 
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INFLUENCE  OF  LIQUEFACTION  ON  LASER  ABLATION: 
DRILLING  DEPTH  AND  TARGET  RECOIL 


A.D.  Zweig 

Wellman  Laboratories  of  Photomedicine 
Massachusetts  General  Hospital 
Boston,  MA  02114 

ABSTRACT 

Ablation  of  materials  by  highly-absorbed  laser  pulses  (duration  200  ps)  is  modelled 
using  a  steady-state  approach  that  takes  into  account  target  hquefaction.  Material  is 
removed  by  a  combination  of  evaporation  and  ejection  of  liquid  that  is  caused  by  radial 
evaporation-induced  pressure  gradients  at  the  bottom  of  the  ablation  crater.  As  a  result 
the  evaporation  rate  is  reduced  and  the  drilling  efficiency  enhanced.  The  model  explains 
experimental  drilling-depth  and  recoil-momentum  data 

SUMMARY 

Drilling  depth  and  target  recoil  provide  important  and  easily  measurable  information 
on  laser  ablation.  They  vary  with  the  distance  between  the  surfaw:e  of  the  target  and 
the  waist  of  the  laser  beam.  For  laser  pulses  of  several  100  ps  duration  the  observed 
variation  is  very  similar  for  biological  targets  and  metals,  suggesting  that  the  same 
physical  mechanisms  are  involved  in  ablating  these  materials'’*. 

Targets  illuminated  with  sufficiently  intense,  highly-absorbed  laser  radiation  are 
heated  at  the  surface  and  start  to  evaporate  after  a  ^ort  fraction  of  the  ablating  laser 
pulse.  Because  of  momentum  conservation  the  resulting  violent  vaporization  of  material 
increases  the  pressure  locally.  To  a  first  approximation  the  pressure  is  proportional  to 
inadiance.  The  induced  pressures  vary  spatially  with  the  beam  profile  of  the  ablating 
laser.  Thus  if  the  target  material  liquefies  prior  to  evaporation,  the  pressure-induced 
shear  forces  at  the  surface  push  mederial  in  radial  direction  towards  the  crater  wall. 
There  the  liquid  flow  changes  direction  and  continues  axially,  along  the  wall  and  out 
of  the  hole.  This  mechanical  ejection  mechanism  becomes  more  efficient  as  more  liquid 
moves  in  the  radial  direction.  CSorrespondingly,  drilling  is  enh^mced  by  a  low  viscosity 
of  the  molten  material,  a  large  thickness  of  the  liquid  surface  layer,  and  a  small  spot 
size  of  the  drilling  beam. 

For  a  given  laser  pulse  the  total  drilhng  depth  changes  with  the  distance  from  the 
surface  of  the  target  to  the  waist  of  the  beam.  Under  tight  focusing  conditions  the  me¬ 
chanical  drilling  contribution  becomes  important,  as  can  be  seen  from  the  experimental 
results  presented  in  Fig.  lau  The  solid  cinve  is  the  result  of  a  model  calculation  based  on 
a  steady-state  description  of  drilling^.  The  calculation  took  into  account  the  variation 
of  irradiance  with  distance  from  the  beam  waist  and  used  parameters  corresponding 
to  the  experimental  data.  Also  shown  in  Fig.  la  is  the  calculated  depth  for  purely 
evaporative  drilling  into  a  water  target  (dashed  fine).  As  is  made  clear  by  the  figure. 
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Figure  1:  Influence  of  the  initial  position  of  the  target  surface  on  the  drilling 
mechanism.  Positive  (negative)  positions  correspond  to  a  beam  waist  lo¬ 
cated  outade  (inside)  the  target.  Experimental  data  for  gelatin  containing 
0.84  g/cm^  water  is  indicated  by  filled  circles  (•).  The  short-dashed  lines 
are  to  guide  the  eye.  All  the  holes  were  drilled  with  single  laser  pulses 
(A  =  2.94  /im,  T  =  200  fis,  tJo  =  72  fim.  E,r.  =  40  mJ).  Data  from’, 
(a);  Drilling  depth  versus  waist-target  distance.  The  sohd  line  is  the  result 
of  a  model  calculation  that  takes  into  account  evaporation  and  ejection;  the 
dashed  line  shows  the  calculated  depth  for  a  purely  evaporative  drilling,  (b): 
Coupling  coefficient,  C,  versus  waist-target  distance.  C  is  equal  to  the  total 
recml  momentum  normalized  by  the  laser  pulse  energy.  Note  that  C  has  a 
relative  mifiimum  when  the  drilling  depth  has  a  maximum. 


the  experimental  values  are  too  large  to  be  caused  by  evaporation  only.  The  enhanced 
drilling  depth  originates  frxim  the  ejection  mechanism. 

The  reo(^  momentum,  caused  by  the  ejected  material  during  an  ablation  process  on 
a  target,  provides  another  simple  means  of  characterizing  the  interaction  experimentally. 
The  total  recoil  momentiun  divided  by 'the  pulse  energy  is  called  the  mechanical  coupling 
ooeffideiit,  C,  and  has  units  of  dyn  ■  s/J.  In  Fig.  lb  its  dependence  on  the  beam  waist 
to  target  distance  is  shown.  The  appearance  of  that  curve  is  similar  for  a  variety  of 
biological  targets  as  well  as  metals.  It  can  be  explained  by  the  generation  of  liquid 
material  during  ablation. 

N^ecting  the  contribution  of  liquid  ejecta  to  the  tot2Ll  target  recoil,  the  coupling 
coefficient  can  be  estimated  roughly  by  assuming  that  all  the  evapioiated  mateiial  leaves 
the  target  at  the  speed  of  sound.  This  picture  is  valid  only  if  the  induced  gas  flow  is 
one-  dimen sirmal  Fcr  high  aspect-ratio  holes  this  can  be  expected  despite  the  long  pulse 
durations  used  (200  fts). 

Therefore  the  recoil  curves  can  be  interpreted  by  considering  the  influence  of  aspect- 
ratio  and  total  mass  evaporated  as  follows:  For  shallow  boles  the  vapor  flows  mainly 
parallel  to  the  target  surface,  with  but  a  small  velocity  component  perpendicular  to  the 
surface  of  the  target.  As  the  radial  momentum  components  cancel  each  other  because 


